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Algorithm  for  the  Calculation  of 
Absorption  Coofficient-Prosiure 
Broodonod  Moiv^culor  Transitions 


I.  INTRODdCTIUN 

ThU  raport  d«aorlbM  an  algorithm  and  aaaoctated  program  to  perform  the 
convolution  of  a Una  ehape  function  with  an  array  of  abeorptlon  tranaltlona  for 
which  the  wavenumber  valuea,  the  inteneltlee,  and  the  halfwldthe  are  apeotflod. 

The  algorithm  haa  been  developed  with  the  objectlvee  of  minimizing  the  number  of 
operatlone  to  perform  the  convolution  and  of  developing  a program  having  reaaon* 
able  computer  etorage  requlremente.  The  decreaue  In  computational  operatlone 
achieved  by  the  preaent  method  may  be  utilized  to  decrenae  computer  time  require - 
menta  for  current  problemu  or,  If  neceusary,  to  perform  more  extenatve  calcula" 
tlona  that  ware  not  previously  considered  feaalble. 

The  program  described  In  thla  report  uses  a given  line  shape  lUnotton  (the 
Lorentf.  hmctlon)  which  la  assumed  to  be  Independent  of  molecular  species,  The 
algorithm  enablas  the  calculation  of  the  abaorptlon  coefficient  at  equally  spaced 
output  wavenumber  Increments  over  a epeolfled  wavenumber  region,  A uniform 
path  la  aaaumed  for  which  the  tamperature,  preaaure,  and  absorber  amounta  are 
specified.  The  program  haa  been  written  to  be  compatible  with  the  AFQL  line 
Hating,  and  performs  the  function  of  program  UN  prevloualy  utilized  In  HITRAN 
(Received  for  publication  3 1 July  1977) 

1.  MoClatchay,  R.A,,  Benedict,  19.3.,  Clough,  S.A,,  Burch,  D.G,,  Calfet, 

R.F. , Foa,  K,.  Rothman,  L.  S. . and  Oaring,  J,  S,  (1973)  AFCRl,  Atmos  - 
pherlo  Abaorptlon  Line  Parameters  Comellatton.  Rsport  No.  afcrl-TR- 
ghifeonm^^^^  Pipers."  Wo.  >134. 
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calciilutlons.  The  algorithm  ban  been  developed  with  the  goal  of  keeping  the  error 
in  the  HbHorptlon  coefficient  to  within  0.  1 percent. 

The  concept  that  has  been  developed  involvcH  the  deoompo.sition  of  the  given 
line  shape  function  into  ihroe  functions  which  are  independently  convolved  with  the 
parameters  for  each  spectral  lino.  The  flr.st  function  is  niirrow  and  describes  tho 
central  portion  of  the  Lorcnt/.lan  line  shape,  the  second  function  is  broader  describ- 
ing tho  Intermodlute  region  of  the  shape,  and  tho  third  function  Is  broad  and  may 
bo  considered  to  describe  the  "near'*  wing  contribution  of  tho  line.  The  functions 
have  been  carefully  cho.son  so  that  the  required  spectral  Information  is  retained  as 
each  function  Is  sampled  at  the  same  number  of  points  by  the  convolution  process. 
The  spectrum  Is  reconstructed  at  tho  final  stage  by  properly  Interpolating  tho  three 
ftmotlons  Into  the  output  array  of  absorption  coefficients.  Also  included  in  this 
report  is  a program  that  performs  a subsequent  Interpolation  providing  nn  essen- 
tially continuous  resultant  spectrum  without  requiring  the  convolution  calculations 
to  be  performed  at  more  output  points  than  are  required. 

2.  DDCOMPOHITKIN  OK  TIIK  I.OHKM'/  MM;  SIIAPK 

The  Lorentz  shape  is  generally  appltoable  for  pressure  broadened  lines  for 
which  the  transition  wavenumber  value  Is  large  compared  to  the  lino  halfwldth,  and 
has  been  dlscuased  by  many  authora.  ^ The  abaorptlon  ooefflcient,  us  a 

function  of  wavenumber,  v,  resulting  from  the  convolution  of  the  Lorentz  line  shape 
with  a molecular  Iranaltlon  of  wavenumber  Inlenalty  Sj,  and  width  (halfwldth 
at  half  maximum  I HWHM)  la  given  by 


and  the  total  abaorptlon  coefficient  la  given  by 


Mv)  • ^ A^jiv)  . (2) 
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For  computational  purpoEie.s,  it  becomes  advantageous  to  define  a dimenaionlcsB 
variable  z which  for  a single  line  is  defined  as 


where  r.  is  the  wavenumber  difference  from  the  line  center,  In  terms  of  the 
halfwldth  of  the  tranaition.  The  absorption  coefficient  for  a single  line  In  terms 
of  z is 


“ 1 7^7  ‘ (I)  (4) 

where  the  Lorentz  function  Is 

Uzl  ■ . (5> 

1 + 

This  definition  has  the  computational  advantage  that  L(z)  is  cxpllcltLy  Independent 
of  halfwldth  and  may  be  tabulated  as  a numerical  function  with  argument  z. 

The  function  L<z)  has  pathologically  slow  convergence  to  zero  as  has  been  dls- 
cuaaed  by  Qoody^  and  others.  In  order  to  avoid  computation  of  the  Lorentz  function 
at  values  of  the  argument  for  which  the  function  is  slowly  varying,  the  function  has 
been  decomposed  into  three  bounded  regions  (0  £ | z|  < Z^,  f > 1,  3).  We  have 
chosen  • i,  Z,^  • 12,  and  Z,^  = 48  resulting  In  three  functions  with  similar  func- 
tional behavior  over  the  bounds  0 to  3 halfwidths,  0 to  12  halfwidths  and  0 to  48 
halfwidths  respectively. 

The  decomposition  has  been  performed  in  a progressive  manner  by  creating 
the  first  function  (denoted  as  the  "fast  function, " (XF(z)>  as  the  difference  between 
the  Lorentz  function  and  cn  even  quartlc  function  Q^<z)  chosen  such  that  the  value, 
the  first  derivative  and  the  second  derivative  of  XP(z)  are  zero  at  the  boundary 
Zy  XF(z)  is  given  by 

XF(z)  ■ l/»  (L'(z) -Qj(i))  for  0s|z!a:Z^  (6) 

where 

L'(z)  • tL(z)  ■ — ^ (7) 

I + z' 


1 


(8) 


I 


i 


Q,<z)  = aj  + bi/2  4 


and 


Cj-l/d  + Z^)^  , bj  ■ -C|(l  + 3 Z^)  and  a^  ■ Cj(l  + 3Z*  + 3Z^) 


For  x Z|  * 3 wa  obtain 

°l“~V  • ‘*1""“^  ' ‘ 

‘ 10^  ' 10^  ^ 10“* 


(9) 


(10) 


Conaequently,  wa  have  for  tho  "faat  funotlon" 
XF(*)  ■ l/»r  <1/(1  + z*)  - (ij  + bj  z^) 

XF(«)  ■ 0 


0 S 1 zl  *3  , 

3 S 1 *1  . 


(ID 

(12) 


Thfl  aacond  funotlon,  the  "alow  funotlon,"  XS(z),  la  oonatructed  In  a almllar 
mannar  by  defining  an  even  quartlc  (^notion,  domain  3 s |z|  ss  12 

whloh  matchea  tho  Lorontzlan  at  e ■ ■ Z^  ■ 12.  We  obtain  by  performing  opera- 

tlona  almllar  to  thoaa  above 


X3(b)  • l/»  (L'(z)  - CJj(z))  for  3 s | z]  a 12  (13) 

with  the  Goafflclenta  for  Q2(>)  given  by 


Oj  ' 74^3  • ‘’2  ‘ • ®3  • **“*  *2 

« 63641  • Oj  . 

(14) 

The  "alow  function"  la  than  defined  ae 

XS(b)  - l/»  (Qj(b)  -Qj(a)) 

0 as  Izl  a 3 , 

(15) 

X3(b)  ■ 1/v  (1/(1  + a*)  - (6j,  + bjB*  + Og  z*)) 

3 a;  |x|  as  13  , 

(16) 

XS(B)  - 0 

12  ii  |z|  , 

(17) 

The  third  funotlon,  the  "vary  alow  Amotion, " XVS(z),  la  defined  aa  the 
Lorentalan  for  13  rt  U|  ^ 4B  and  aa  Q^^a)  for  0 « |z|  as  13, 
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The  following  tabular  form  of  the  decomposHlon  Into  the  respective  domains 


will  moke  the 

proce.SH 

more  clean 

FUNCTION 

0 s |zl 

s 3 

3 £ |z|  ^ 12 

12  s |z|  48 

XF(z) 

i (L'(Z) 

- g^(z)) 

0 

0 

XS(z) 

i (Ql(z) 

-Q2(z)) 

^ (L'(z)  - Qjfz)) 

0 

XVSlz) 

i Q2(*) 

i a,(.) 

Note  that  the  three  functions  sum  to  the  Lorentzlun  in  each  domain  and  that  the 
functions  are  continuous  In  value,  first  and  second  derivative  across  the  boundaries 
of  the  domains. 

While  the  choices  of  the  functions  and  the  boundaries  may  appear  to  be  some- 
what arbitrary,  an  examination  of  the  result  demonstrates  the  advantage  of  the 
Imposed  conditional  Figure  1 Is  a plot  of  the  three  functions  from  -4B  to  4-4B  half- 
widths and  Indicates  the  dependence  of  XP(z),  XS(z),  and  XVS(z)  as  a function  of 
the  variable  z.  It  should  be  noted  that  40,71  percent  of  the  integral  of  the  Lorentz 
function  la  contained  In  XF,  42.47  percent  in  XS,  and  16.  50  percent  In  XVS.  The 
remaining  1.32  percent  Is  outside  the  48  halfwldth  cutoff. 

The  three  functions  are  represented  in  Figure  2 with  the  horizontal  axis  deter- 
mined by  the  boundary  of  each  function  and  the  vertical  axis  chosen  appropriately 
for  each  function.  The  three  functions  in  thlu  representation  are  seen  to  be  func- 
tionally similar  and  as  a consequence,  it  is  possible  to  sample  each  function  at  the 
same  sampling  Interval  to  perform  the  convolution.  One  hundred  and  fifty  one  (151) 
values  of  each  symmetric  function  are  stored  In  the  computer  program,  and  It  Is 
these  values  that  are  plotted  In  Figure  2 as  Indicated  by  the  Index  on  the  top  hori- 
zontal axis.  The  tabulation  of  thla  number  of  values  for  the  functions  precludes 
the  necessity  for  Interpolation  tn  the  convolution  calculation.  The  valves  of  the 
functions  are  tabulated  in  Table  1. 
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Figure  2.  Functional  Behavior  of  the  Decotnpoaed  Functloni 
Over  Their  Renpective  Limits.  Plots  are  direct  plots  of  the 
functions  in  terms  of  the  storage  index  of  the  pro^am.  The 
fast  function  (XF)  is  on  the  top,  the  slow  function  (XS)  is  in 
the  middle,  and  the  very  slow  functions  (XVS)  are  on  the 
bottom 
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Table  I,  Tabulation  of  the  Decomposed  Functions  Used  to  Describe  the  Lorentz 


Function.  The  variable  z Is  defined  as  the  wavenumber  difference  from  the  line 
center  In  terms  of  the  halfwldth,  a,  that  Is,  z = l<r  - v )/ol.  Plots  of  the  three 
tabulated  functions  appea r In  Figure  2 and  Figure  3 ° 
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.8842-03 

SO 

.880 

1.8497-01 

2.320 

6. 

1212-02 

9.2  80 

.3 

•4211-03 

31 

.6  00 

1.9107-01 

2.490 

3 • 

9202-02 

9.600 

3 

.2612-03 

32 

.6  20 

t. 4707-01 

2.460 

3. 

722E-02 

9.920 

3 

.1032-03 

33 

.6  40 

1.4327-01 

2.860 

3. 

9272-02 

10.24l< 

2 

.9482-03 

34 

. 660 

1.  3937-01 

2.640 

3. 

338E-02 

10.864 

2 

.7972-03 

38 

.660 

1.3887-01 

2.720 

3. 

183E-02 

10.880 

2 

.6822-03 

36 

.700 

1.3177-01 

2.800 

2. 

976E-02 

11.200 

2 

.8122-03 

37 

.7  20 

1.2797-01 

2.880 

2. 

806E-02 

11.820 

2 

. 3792-03 

36 

.740 

1.  2427-01 

2.960 

2. 

6462-02 

11.840 

2 

.2882-03 

38 

.7  60 

1.2087-01 

3.040 

2. 

4982-02 

12.160 

2 

.1302-03 

4? 

.7  80 

1.1697-01 

3.120 

2. 

3842-02 

12.400 

2 

.0312-03 

41 

.880 

1.  1347-Dl 

3,294 

2. 

223fc-02 

12.800 

1 

.9312-03 

42 

.8  20 

1.0997-01 

9.280 

2. 

lOOE-02 

13.120 

1 

.8392-03 

43 

. 840 

1,0687-01 

3.360 

1. 

9*62-02 

13.440 

1 

. 7522-03 

44 

.860 

1.0317-01 

3.440 

!• 

878E-02 

13.760 

1 

.6722-03 

4^ 

.880 

9.9847-02 

3.870 

1. 

7782-02 

14,080 

1 

.8982-03 

46 

. 900 

9.6617-02 

3.600 

1. 

6832-02 

14.400 

1 

.82*2-03 

47 

.920 

9.  3497-02 

3.660 

1. 

8942-02 

14,720 

1 

.4622-03 

48 

.940 

9.0387-02 

3.768 

1. 

81U-02 

18.040 

1 

.4012-43 

49 

. 960 

8,7337-02 

3.640 

1« 

4322-02 

18.360 

1 

.3432-03 

80 

.980 

0.43T7-02 

3.97s 

1. 

3882-02 

18.680 

1 

.2*92-03 

81 

1.0  90 

6.199C-02 

4.900 

1. 

2886-02 

18.000 

1 

• 2392-03 

82 

1.0  20 

7. 8677-02 

4.080 

1. 

2222-02 

16.320 

1 

.1912-03 

83 

1.040 

7.  9927-02 

4.160 

1. 

180E-02 

16.660 

1 

.1482-03 

84 

1.080 

7.1247-02 

4.240 

1. 

1012-02 

16.960 

1 

.1032-03 

88 

1.0  80 

7.0617-02 

4.320 

1. 

0466-02 

17.280 

1 

.0622-03 

86 

1.100 

6.8097-02 

4.439 

9. 

9306-03 

17,600 

1 

.0242-03 

87 

1.120 

6.8617-02 

4.680 

9. 

632E-03 

17.920 

9 

.8822-04 

88 

1.140 

8.3207-02 

4.868 

8. 

9606-03 

16.260 

9 

.8392-04 

84 

1.184 

8.0867-02 

4.640 

8, 

9132-03 

18.960 

9 

.2142-04 

60 

1.180 

8.0887-02 

6.720 

8. 

0896-03 

ll.liO 

8 

.9012-04 
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Table  1,  Tabulation  of  the  Decomposed  Functions  Used  to  Describe  the  Lorentz 


Function.  The  variable  a is  defined  as  the  wavenumber  difference  from  the  line 
center  in  terms  of  tlie  halfwldth,  a,  that  is,  7.  = liv  - Plots  of  the  three 

tabulated  functions  appear  in  Figure  2 and  Figure  3 (Cont. ) 


INDEX 

z 

••FAST’' 
rUFJCT  ION 

FN 

z 

"St OH" 
FUNCTION 

FN 

7 

‘•vert  SLOW" 
FUNTTION 

F N 

61 

1.200 

5.6372-02 

6.600 

7.6672-03 

19.200 

6 .6112-06 

6? 

1.2  20 

5.6222-02 

6.660 

7,3092-03 

19. 520 

6.3322-06 

63 

1.260 

5.21  3E-02 

6.960 

6.9622-03 

19,  H40 

5. 0662-06 

64 

t.?60 

5.01 02-02 

5.  060 

6.5972-03 

20.160 

7 ,3132-06 

65 

1.2t0 

6.6132-02 

5.120 

6.2692-03 

20.660 

7 .5712-06 

66 

1.300 

6.6222-02 

5.200 

5.9502-03 

20.600 

67 

1.320 

6.637E-02 

5.260 

5.6612-03 

21.120 

7.12UE-06 

66 

1.360 

6,2582-02 

5.360 

5.3792-03 

21.660 

6 .9102-06 

69 

1.  360 

6.0662-02 

5.660 

5.1112-03 

21.760 

6.  7062-06 

70 

1.3B0 

3.915E-a2 

5.520 

6.0582-03 

22.060 

6 .5162-06 

7l 

1.600 

3.7522-02 

. .600 

4.6112-03 

22.600 

6 .3312-06 

1.620 

3.5962-02 

5.630 

4.3792-03 

22.720 

6.1552-06 

73 

1.660 

3.6612-02 

5.760 

4.1582-03 

23.060 

9.9352-06 

74 

1.660 

S.293r-o» 

5.860 

3.9672-0  3 

23. 360 

5 .8322-06 

75 

1.660 

3.1502-02 

5.920 

3.7462-03 

23.680 

666E>04 

76 

t.too 

3.0122-02 

6,000 

3.5552-03 

76,000 

5 .5172-06 

77 

1.6  20 

2.0792-02 

6.080 

3.  3722-0  3 

26. 320 

5.3732-06 

76 

1.  660 

2.76  92-02 

6.168 

3.1982-03 

26.660 

5. 2362-06 

79 

1.660 

2.6252-02 

6.260 

3.0322-03 

26.960 

5.1012-06 

60 

i.tto 

2.5062-02 

6.120 

2.8732-03 

25.260 

6 .97  3 2-06 

61 

1.  600 

2,3882-02 

6.600 

7.7222-03 

25.  600 

6.  6502-04 

62 

1. 620 

2,2762-02 

6.6  00 

2.5782-03 

25.920 

6.7312-06 

63 

1.660 

2.1682-02 

6.960 

2.6612-03 

26,260 

6 .6162-06 

64 

1.660 

2.0662-02 

6.660 

2.3692-03 

26. 560 

6 .5062-06 

65 

1.6t0 

1.96  32-02 

6.720 

2.1862-03 

26.660 

6 .3992-06 

66 

1.700 

1.8662-02 

6.808 

2.0652-03 

27.200 

6.2972-06 

67 

1.720 

1.7732-02 

6.880 

1.9512-03 

27.  920 

6,1972-06 

66 

1.  760 

1.6662-02 

6.960 

1.8632-03 

27.660 

6.1022-04 

69 

1.760 

1.3972-02 

7.060 

1.7S92>03 

26.160 

4.0092-06 

90 

1 « 760 

1.5162-02 

7.120 

1.6612-03 

26.660 

3 .9202-06 

91 

1.600 

1,6352-02 

7.200 

1.5672-03 

26,  800 

3.3332-06 

9? 

1.620 

1.35  82-02 

7.208 

t.6J7E-03 

29.120 

3 ,7692-06 

93 

t.evc 

1.2862-02 

7.360 

1.3722-03 

29.660 

94 

1.660 

1.2162-02 

7,660 

1.2912-03 

29.  760 

3*  99jl,'*04 

95 

1.660 

1.1662-02 

7.520 

1.2132-03 

30.060 

3 .5162-06 

96 

1.900 

1.0612-02 

7.  600 

1.1602-03 

30.400 

3 .wi,12-06 

97 

1.920 

1.0192-02 

7,680 

1.0702-03 

30  . 720 

3 •349E-04 

96 

1.  960 

9,5962-03 

7.760 

1.0032-03 

31. 660 

3. 3002-06 

99 

1«960 

9.0262-03 

7,860 

9.3972-06 

31.360 

3 .2332-06 

lao 

1.960 

6.6792-03 

7.920 

8.7962-06 

31.680 

3 .1652-04 

101 

2.0  00 

7.95  82  - 03 

8.030 

8.2262-04 

32. COO 

3»105E-I)4 

102 

7.020 

7.66  02*03 

0.0  80 

7.6022-06 

32.320 

3 •044E-04 

103 

2.  060 

6.9852-03 

8.160 

7.1672-06 

32.660 

104 

2.0  60 

6.5322-03 

6.260 

6.6792-06 

32.960 

?a9Z7C-04 

105 

2.060 

6.1012-03 

0.321! 

6,2172-06 

33  .260 

z •eMC-04 

106 

2. 100 

5,6902-03 

8.690 

9.7792-06 

33. 600 

2tBlFE-04 

lor 

2.  120 

5,2992-03 

6.680 

9.3652-06 

33.920 

7B4E-04 

100 

2t  140 

6.9282-03 

6.560 

6.9732-06 

36,260 

2 f7l3E-04 

loa 

2.160 

6.5752-03 

8.660 

4.  60  32-04 

36,560 

Z •663E»a4 

110 

2.  ISO 

6.2612-03 

a .7  20 

4. >942-06 

36. 880 

2«414E«04 

111 

2.200 

3.9262-03 

8.  80  0 

3.9292-06 

35.200 

E«56?C-04 

112 

2.220 

3.8262-03 

8.880 

3,  6192-0  6 

35.520 

2 .5212-06 

113 

2.260 

3jl6ie>D3 

8*960 

3.  3232-06 

35.  660 

2.6762-06 

114 

2.260 

3 . 0 7 32  - 03 

9.060 

3.0692-96 

36.160 

2.6332-06 

in 

2.260 

2.8212-03 

9*128 

2.  7912-04 

,36.680 

2 .3902-06 

116 

2.300 

2.90  62-03 

9.208 

2.9902-04 

36,600 

2«349E*04 

117 

Z«  370 

2.3812-03 

9.280 

2.3242-04 

37,120 

2.3082-06 

lit 

2.360 

2.1922-03 

9.360 

2,1132-06 

37.660 

2 .2092-06 

in 

2.360 

1.9982-03 

9,668 

1.9162-06 

37,  760 

. .2312-06 

120 

7 • 3 AO 

1.77  32-03 

9.920 

1.7322-06 

38.010 

2.1942-06 
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Table  1,  Tabulation  of  the  Decomposed  Functions  Used  to  Describe  the  Lorentz 
Function.  The  variable  z Is  defined  as  the  wavenumber  difference  from  the  line 
center  In  terrn.s  of  the  halfwldth,  a,  that  is,  z [(t'  - Plots  of  the  three 

tabulated  functions  appear  in  Ftjpiro  2 and  Figure  3 (Coni.) 


INOEK 

7 

"7*57“ 

FUNCTION 

FN 

2 

“SLOW" 

FUNCTION 

FN 

: 

“VERY  StOM" 
FUNCTION 
FN 

IJl 

2.400 

1.6D2E-03 

9.600 

1.561E-04 

38.400 

2.1572-04 

12? 

2.420 

1.442F-01 

9.600 

1.40SE-04 

36.720 

2.1222-04 

1?] 

2.44(1 

1.294E-03 

9.760 

1.256E-a4 

39.040 

2.007E-04 

1?4 

2.460 

1.157e-03 

9.940 

1.121E-0'! 

39.360 

2.Q5JE-04 

125 

2.400 

1.0300-03 

9.920 

9.955E-05 

39.600 

2.b2UE-04 

126 

2.500 

9.131E-04 

10.000 

I.006E-05 

40.000 

1 .90SE-04 

12? 

2.520 

2.  055E-fl4 

10.090 

7.753E*0S 

40. 320 

1.957E-04 

121 

2.940 

7.070E-04 

10. 160 

6.791E-05 

40.640 

1 .926E-04 

125 

2.560 

6.170E-04 

10.240 

5.915E-U5 

40.960 

1.996E-04 

139 

2.5  00 

5.392E-04 

10,320 

5.121E-05 

41.20r 

1.9671-04 

131 

2.600 

4.610F-04 

10.400 

4.403E-05 

41.601 

1 .039E-04 

132 

2.620 

3.942E-n4 

10.400 

3.  750E-05 

41.920 

1.81JE-04 

133 

2.640 

2.343E'J4 

10.560 

3.102E-05 

42.240 

1.799E-04 

134 

2,660 

2.909E-C‘- 

10.640 

2.6b9E-05 

42.560 

1 .7561-04 

135 

2,600 

2.336E-04 

10.720 

2.  215E-D5 

42.000 

1 .7I0E-04 

136 

2.730 

1.<)20E-04 

10.000 

1.0197-05 

43.200 

1.  70  5E-04 

137 

2.7  20 

1.0070-04 

10.000 

1.472E-05 

43.520 

1.6O9E-04 

13t 

2.740 

1.244E-04 

10.960 

1.174C-05 

43.  040 

1 .655E-04 

135 

2.  7 60 

9.70  9E-05 

11.940 

9. 1942-06 

44.160 

1.631E-04 

140 

2.700 

7.499E-00 

11.120 

7.054E-D6 

44.400 

1.6091-04 

141 

2.9  00 

O.620E*bO 

11.200 

5.2SaE-06 

44.000 

1.596E-04 

14? 

2.  020 

4.089E-00 

11.260 

3.034E-06 

45. 120 

1.563E-04 

143 

2.9  40 

2.  S60E-00 

11.360 

2.699E-06 

45.440 

1.941E-04 

144 

2.960 

1.915E-05 

11.440 

1.791E-06 

45.76(1 

1.519E-D4 

146 

2.9  90 

1.2C3f-fl5 

11.520 

1.124C-00 

46.090 

1.499E-04 

146 

2.900 

6.947^-06 

11.600 

6.4792-07 

46.400 

1.4797-04 

147 

2.9  20 

3.049f-06 

11.600 

3.  306E-U7 

46.720 

1.459E-D4 

140 

2.940 

1.494£-06 

11.76(1 

1.  590E-D7 

47.  040 

1.  430E-(i4 

149 

7.960 

4.419E-07 

11.040 

4.104E-00 

47.360 

1.419E-04 

190 

2.990 

5.011E-P0 

11.920 

5.113E-09 

47.690 

1.40JE-04 

151 

3.000 

C. 

12.900 

D. 

49. ODD 

1.391E-D4 

The  ratlunallsjatlon  for  terminating  the  convolution  at  48  halfwidths  is  based 
on  three  considerations!  (1)  evidence  that  the  Iliie  shape  deviates  from  Lorentz 
behavior  for  (r  - v ) of  the  order  of  6 cm"',  correspondlnr  to  i » 46  for  a typical 

O r fl  Q 

atmoapheric  line  at  atmospheric  pressure  (Winters  et  al:~'  Holstein;  Breene; 
Burch;^)  (2)  the  value  of  the  Lorentzian  at  48  halfwidths  is  1. 38  X lO"^  which  is 


5.  Winters,  13.  H.,  Silverman,  S.,  and  Banedict,  W.  S.  (1964)  Line  shape  in  the 

wing  beyond  the  band  head  of  the  4. 3 p band  of  CO,,  J,  Quant.  Spect,  Rad. 
Trans,  4iB27,  * 

6.  Holstein,  T,  (1660)  Prassure  broadanlng  of  spectral  lines,  Phys,  Rev.  79:744; 

sea  also  L.  Spltzer  Phys.  Hev.  Mt348  (.  ^0). 

7.  Burch,  D.E.,  Orynak,  D.A.,  Patty,  R.R.  and  Bart ky,  C.E.  (1608)  Ibft. 

Shapes  of  Colltalon  Broadened  CO„  Lines,  Phllco  Ford  Corp..  Aeronutronlo 
Heport  t>3a03.  * 
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Muri'lclently  nmnll  for  nioHt  calciiLatlon.s,  and  (3)  calculation  beyond  46  hairwldth.s 
roqulrea  an  unwarranted  Increase  In  I'omputatlonal  effort.  Although  the  vnhio  of 
the  ^oroiitv,  function  may  justify  termination  of  the  calculation  at  46  halfwldths,  It 

must  not  be  forgotten  (hat  1.32  porennf  of  the  Integrated  absorption  Is  1,'eyond  this 

2 

region.  Thin  Is  a result  .o'  iho  previously  mentioned  slow  convergence  of  I// 

2 fl 

jQoody,  and  Townes  and  Srhawlow  ).  It  Is  suggested  that  calculations  of  the  far 
wings  of  the  lines  can  inure  appropriately  be  done  in  a parametrised  tabulation  of 
what  Is  currently  referred  to  as  the  continuum  for  each  of  the  molecular  species 
of  Interest.  Thu  contribution  of  this  continuum  function  to  the  absorption  coeffi- 
cient may  be  calculated  at  large  sampling  intervals  In  the  wavenumber  domain  with 
negligible  effect  on  computational  time. 

3.  SAMPIJNli  INTKHVAI. 


General  discussions  of  the  proper  choice  of  the  sampling  Interval  appear  In 
the  literature  for  the  type  of  "line  by  line"  calculatlone  deacrlbed  In  this 
report.  ^ These  references  are  primarily  concerned  with  the  effect  of 
the  sampling  Interval  on  the  integral  of  the  absorption  coefficient  or  ebaorptance 
over  a specified  wavenumber  Interval.  That  approach  Is  particularly  appropriate 
to  the  case  In  which  an  Instrumental  scanning  function  that  Is  broad  compared  to 
‘.he  spectral  line  widths,  is  being  convolved  with  the  true  line  spectrum.  The 
present  discussion  is  concerned  with  the  development  of  a quantitative  evaluation 
of  the  error  introduced  In  the  reconstruction  of  the  spectrum  In  the  Interval  be- 
tween sampled  points. 

An  approach  that  provides  direct  Insight  Into  the  problem  of  determining  the 
sampling  Interval  is  that  obtained  from  Information  theory  and  the  Nyqulat  sampling 
theorem.  Simple  and  readily  comprehensible  discussions  of  this  ere  provided  by 


6.  Townes,  C,  H.,  and  Sohawlow.  A.  L.  (1956)  Microwave  Spectroscopy.  McOrew 
HUl-BookCo. , New  York.  

9.  Orayson,  D.R.  (1967)  The  Calculation  of  Long  Wave  Radiative  Transfer  In 
Planetary  Atmoepheree,  Report  No.  07664^ -T.  UnlvereUy  b/ Mlchlgim. 

Ann  Arbor. 

10.  Kyle,  T.Q,  (1968)  Net  interval  for  calculating  abaorptlon  epeotra.  J,  Opt. 

Soc,  Am.  681192.  ■* ^ 

11.  Soott,  N.  A,  (1974)  A direct  method  of  computation  of  the  tranamlsslon  func- 

tion of  an  inhomogeneoue  gaaeoua  medium,  Ii  Description  of  the  method, 

J,  Quant.  Spect.  Rad,  Trana.  141681. 

18.  Kunda.  V.Q.,  and  MoOuira,  W. C.  (1974)  Dlract  integration  transmlttonca 
modal,  J.  Quant,  Spect,  Red,  Trena.  IdiBOS. 
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13  14 

Brl((ham  and  Bracewell.  The  »ampUng  theorem  atatex  that  it  the  Kourler 
tranaform  of  a function  la  zero  for  all  frequcnctex  greater  than  frociuency  then 
the  rontinuourt  function  g(x)  can  bo  uniquely  determined  from  ita  sampled  valuca, 
(nAx),  if  the  sampling  Interval  la  chosen  as 

•-■d-;-  O'" 


In  particular,  the  contlnuoua  function  g<x)  can  be  reconstructed  for  all  x by  con- 
volving the  contlnuoua  Interpolating  function 


with  the  sampled  function  g^  (nAx)  where  n la  Integer.  That  la, 

g(x)  • g^  <nax)  * a(x)  (21) 

or  explicitly  for  the  discrete  convolution  we  have 


to 

g(x)  • ^ g^  (MX) 
n«-« 


(23) 


In  the  present  case,  the  Fourier  traneform  of  the  Lorents  function  does  not 
have  a frequency,  f^,  beyond  which  the  tranaform  la  zero.  Consequently  by  sam- 
pling at  dleorete  Intervals,  It  Is  not  possible  to  exactly  reconstruct  the  contlnuoua 
function.  However,  by  properly  choosing  the  Interval,  the  difference  between  the 
exact  and  reconstructed  functions  can  be  held  within  a predetermined  limit.  The 
Lorents  Amotion  Is  given  as 

g(x)  ■ I afl^  A (33) 

a + X 


13.  Brigham.  B.O.  (1B74)  The  Fast  Fourier  Tranaform,  Englewood  Cliffs. 
New  Jersey.  — — — — 

Bracewell,  R.M.  (19SB)  The  Fourier  Tranaform  and  Ita  ApDlloatlons. 
lloaraw-HlU,  New  YorET 
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where  x - i'  - n-  irt  the  halfwldlh,  and  In  the  wavenumber  value  of  the  tran- 
altlon,  The  Fourier  tranaform  of  n(x)  la 


j^(k(x»  ■ ^(f)  ■>  0 


■2if<vf 


<24) 


where  ^ la  the  Fourier  transform  operator.  If  the  value  at  the  tranaform  la 
aaaumed  to  bo  *ero  beyond  a frequency,  f^,  thla  la  equivalent  to  multiplying  the 
complete  tranaform  by  a rectangular  box._  where 


fl  ^ 


We  may  define  the  aampllng  Interval  ae  before  auch  that 


" af: 


The  Fourier  iranaform  olj^^ (f)  In  the  x domain  la  given  by 

j atn  (t^) 
a(x)  ■ ~ — T'-: 


Ax 


(26) 


(20) 


(27) 


which  la  Identical  to  the  acannlng  function  prescribed  by  the  Nyqulst  sampling 
theorem  (Eq.  (20)). 

In  order  to  evaluate  the  error  caused  by  the  discrete  aampllng  Interval  and 
the  associated  truncation  of  the  function  In  the  Fourier  domain,  calculations  have 
been  made  using  the  Interpolating  function,  Eq,  (27),  convolved  with  the  sampled 
Lorenta  function 

g,  (ndx)  « 4 j (28) 

• o*  + (ndx)* 

as  Indicated  In  Eq.  (22).  For  these  convolution  oaloulatlons,  a more  relevant 
definition  of  the  aampllng  Interval,  Ax,  la  given  In  terms  of  the  halfwldth,  o, 

6*  a . (28) 
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'I'hc  romillH  uro  uvnlUBtcd  ,'or  HumpUng  Intorvals  expreaaecl  In  fructlona,  1/N,  of 
(ho  hulfwldth.  For  thla  definition,  wo  have  for  the  sampled  Forenlz  function 


A quantitative  evaluation  of  the  effect  of  the  aampllng  Interval,  la  the  RMS  devia- 
tion of  the  function  Interpolated  at  three  Intermediate  points  per  Interval  from  the 
function  calculated  exactly  at  the  same  Intermediate  polnta.  The  reaulta  of  these 
calculations  are  given  In  Table  2.  These  reaulta  are  a measure  of  the  amount  of 
Information  lost  due  to  the  choice  of  sampling  interval.  Aa  the  sampling  Interval 
becomes  Infinitesimally  small,  the  KMS  deviation  approachea  xero.  It  Is  aeen 
that  for  an  Interval  of  a/i,  that  Is,  N • 4,  the  RMS  deviation  Is  negligibly  small 
and  essentially  all  of  the  Information  la  retained. 


Table  2.  Effect  of  Sampling  Interval  on  the  Reconstruction  of  the  Lorentx  Function 
Using  (sin  x/x)  Interpolation  and  Pour  Point  Lagrange  Interpolation,  Interpolations 
are  calculated  at  three  equally  spaced  points  between  the  sampled  values,  In  all 
oases  (he  maximum  percent  deviation  using  the  four  point  Lagrange  occurs  In  the 
first  sampled  Interval.  Deviations  are  calculated  as  the  difference  between  the 
Interpolated  value  and  the  exact  value  of  the  function 


Interpolation  Method 

Sampling 

Interval 

aln  X 

X 

Pour  Point  Lagrange 

N 

Ax  n 

OX 

rms  deviation 

ims  deviation 

max  paroent 
deviation 

1 

a 

1. 1 X 10"^ 

3.  6 X 10"^ 

-1.  0 

2 

a 

I 

0.067 

0.20 

-0.  01) 

3 

a 

T 

0.  0039 

0,007 

-0.37 

4 

a 

4 

0.0Q063 

0.026 

-0.  IS 

There  Is  another  consideration  that  must  be  taken  Into  account  In  this  situation. 
Although  ths  Information  la  retainad  for  a aampllng  Inttrval  of  e/4,  tha  computational 
effort  to  retrieve  thla  Information  using  the  (aln  x/x)  Interpolation  la  Urge  and  In 
some  casta  may  ba  of  tha  aama  magnitude  aa  that  required  to  avaluate  the  function 
directly  at  tha  Intarmadtata  Intarvals. 
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In  ordor  tn  ri'trteva  HUfflcU'tU  Inforinatlon  wltli  reasonable  computational 
effort  and  stnv  within  i>  spnclflod  error  criterion  for  the  absorption  coofflcUMi' 
oulculatlon,  n four  point  Interpolation  function  has  been  n^iod  Instead  of  tlie  pro- 
hibitively Ion)/  "sin  x/x"  interpolation.  The  outp\it  points  are  assumed  to  be 
equally  spa<'od.  Valuc.s  for  the  wcl(dit!n(<  factors  necessary  to  perform  the  inter- 
polation at  three  equally  spaced  Intormollatr  points  are  given  by  the  following 
three  equations! 


■ rk  * rsf 

' ^ ■ fiu 

(31) 

- ^ Klxl  1 

K(x  ^ 6x)  - ■pij  g(x  + aflx) 

Ci") 

■5^1 

- B(x  - 6x)  + g(x) 

(33) 

The  Interpolation  scheme  that  has  been  employed  in  this  development  is  due  to 
Lagrange  and  la  well  described  by  Acton.  In  Table  U,  the  maximum  percent 
deviation  in  the  Interpolated  values  Is  Indicated.  This  error  always  occurs  within 
one  or  two  Intervals  from  the  line  center.  For  a sampling  Interval  of  (»/4,  the 
error  la  0. 16  percent.  Although  this  error  at  the  line  center  is  slightly  larger 
than  desirable,  the  added  complication  of  using  a more  roflnod  Interpolation  does 
not  presently  seem  warranted.  Furthermore,  the  percent  error  decreases  rapidly 
away  from  the  line  oentor. 

The  four  point  Interpolation  scheme  Is  used  for  Interpolating  the  "very  slow 
array, " end  the  "slow  array"  Into  the  "fast  array."  This  la  required  ulncc  the 
sampling  orltarton  of  o/4  haa  been  used  for  the  calculation  of  the  three  arrays 
where  the  effective  a la  appropriate  to  each  array. 

Similarly,  for  the  preaentatlon  of  the  final  results,  the  fast  array  Is  Inter- 
polated Into  a plotting  array  in  a separate  program  which  represents  the  final 
result.  Figure  3 deplete  the  effect  of  *hls  Interpolation  scheme  as  compared  to  the 
result  obtained  with  linear  Interpolation.  It  should  be  noted  that  the  error  for  a 
complete  spectrum  can  only  be  the  same  or  less  than  that  for  a single  line.  This 
is  readily  concluded  from  the  superposition  theorem. 

One  final  oonalderatlon  that  provides  Information  on  the  aampllng  Interval,  la 
the  Haylalgh  criterion.  For  two  spectral  lines  of  unit  strength  and  equal  half- 
widths, isparated  by  two  halfwidths,  let  us  determlns  ths  aampllng  Interval 
required  to  exactly  reproduce  the  minimum  between  the  two  apeotrol  lines. 

IS.  Acton,  F.  S.  (1B70)  Numerical  Methods  That  Work.  Now  York. 

le.  Bom,  M.,  and  Wolf,  S.  (1985)  Prliootplee  of  Optlge.  Oxford,  Pergamon  Press. 
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3029, S 3030.0  , 3030.5 

mVE  NUMBER  (cm'') 


Flgurt  3.  RtcoiutruoUon  of  th«  Shap«  of  « Slngl*  Line  of 
3lr«nfth  »,  lulfwldth,  a « 0, 1,  Ualng  • Sampling  Intarval  of 
Ax  ■ o/4  With  a Sampling  Potot  Poiltlonad  at  tha  Spaotral 
Lilna  Gantar.  Tha  uppar  ourva  uaaa  llnaar  Intarpolatlon  of 
tha  convolution  output  and  tha  lowar  ourva  uaaa  tha  4 point 
Lagraaga  Intarpolatlon.  To  amphaalaa  tha  Una  oantar,  tha 
region  from  -6  to  +5  holfwldtha  la  daplntad 


ao 


C^>n«ldci‘  ono  Hnu  cpntorod  at.  x - 0 and  the  other  line  centered  at  x • x . Wo 

o 

have  for  the  two  I,iit'cntz  functloiw 


/ I 1 <•' 

7 2 .T 

» ♦ X 


and 


gjlx)  " 7 -K 2 

“ t <x  - %y 


(34) 


(33) 


The  Fourier  traneform  for  thU  pair  ta  given  by 

^ (f)  " e'^*"^ + e'^'"^coH  (2»x^n  . (30) 

The  Kayleigh  criterion  poetulatee  • 2u  giving  for  the  tranaform 

e (l  + ooe  4»<»f)  . (37) 

To  exactly  reproduce  the  minimum  between  the  doublet,  one  period  of  the  ovollla- 
tory  function  la  required,  giving 

4«of  • 2?  (38) 

o 

ao  that  tha  out  off  (raquancy  la 

fo  “ ^ 

and  the  eampllng  Interval  from  the  Nyqulat  theorem,  la 

flx  • • o . (40) 

o 

Thla  reault  Indlcatca  that  the  eampllng  Interval  of  o/d,  prevlouely  determined,  la 
well  within  the  Rayleigh  criterion  for  reaolvlng  a Lorentalan  doublet  aeparated  by 
two  halfwldtha. 
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t.  APPLICATION  OF  THF  Al.liOHITIIM 


In  the  pi'oviouft  Mectiona,  two  buulc  Mtcpa  for  an  Improved  c-oniputatlonal  pro- 
cedure to  perform  the  convolution  of  the  Lorentz  line  ahitpe  with  a not  of  ebaorptlon 
line  date  have  been  deacrlbedi  (1)  the  decomposition  of  the  line  shape  Into  bounded 
domains,  and  (2)  the  criterion  for  determining  the  sampling  Interval,  In  this 
section  we  describe  the  application  of  the  algorithm  to  obtain  the  absorption  coeffi- 
cient as  a function  of  wavenumber  for  a set  of  spectral  line  data  parameters  that 
are  appropriate  for  the  specified  values  of  temperature,  pressure,  and  column 
density.  Three  separate  convolutions  utilizing  the  fast  (XP'),  the  slow  (X,S),  and 
the  very  slow  (XVS)  functions  are  performed  on  the  data  set,  The  results  of  the 
three  convolutions  for  a specific  set  of  data  are  displayed  In  tranamlttance  in 
Figure  4.  For  this  case  the  average  halfwldth  Is  of  the  order  of  0.  04  cm'^  so 
that  the  sampling  interval  used  Is  Ax  » aji  • 0. 01  cm"*.  Each  convolution  function 
Is  sampled  nt  approximately  24  values  corresponding  to  Intervals  of  0,  01  cm'^  for 
the  fast  function,  0.04  cm'^  for  the  slow  function,  and  0. 16  cm"^  for  the  very 
alow  function.  Seventy-two  (72)  points  have  been  sampled  to  reproduce  the 
Lorentzlan  from  -48  to  448  helfwldths  Instead  of  384  (4  X 08)  which  would  be  re- 
quired If  the  function  were  sampled  at  equal  Intervals  of  »/4  (0.01  om'S.  Conse- 
quently, the  algorithm  as  applied  In  this  case  provides  a computational  savings  of 
a factor  b 1/3.  Tha  actual  saving  using  an  optimised  version  of  the  conventional 
method  compared  to  the  present  method  is  s factor  of  10.  The  additional  factor  of 
two  Is  due  to  the  use  of  tabulated  rather  then  calculated  functions  and  to  reduced 
overt,  sad  time  in  setting  up  the  Inner  convolution  loop.  The  results  of  the  two 
computational  methoda  ere  given  In  Figure  5 and  no  differences  are  dlscernablc. 

U should  be  noted  that  the  lower  spectrum  of  Figure  5 Is  equivalent  to  the  product 
of  the  three  functione  preaented  In  Figure  4,  since  the  funotiuna  are  all  on  a 
tranamlttance  scale. 

Some  further  conelderstlon  needs  to  be  given  to  the  method  described  here 
with  respect  to  the  theoretical  computational  gain.  If  the  assumption  Is  made  that 
the  line  shape  function  con  be  decomposed  Into  M funotlons,  and  that  the  sampling 
Interval  Is  (ha  same  for  the  flrit  decomposed  function  as  for  the  total  Lorentz 
funotlon,  a/4,  than  the  gain,  Q,  Is 
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In  thla  expression,  Is  Ihe  total  number  of  halfwidths  over  which  the  function 

iB  to  be  computed,  Ib  the  nunnber  of  halfwidths  apannod  by  the  m'th  decom- 
poaed  function,  and  I*  ratio  of  the  sampling  Interval  of  the  first  func- 

tion to  that  required  for  the  m'th  function.  In  the  present  case  for  -•  48,  we 
have  I 


48 


O'  1 + 12  < 


48 


"h^ 


48 

T 


4 


(42) 


There  are  aome  constraints  that  are  not  explicit  In  the  above  expression;  In  par- 
ticular, there  are  llmitatlona  on  how  small  N„  can  be  and  still  have  the  function 
be  reconstructable  at  the  sampling  interval.  The  algorithm  provides  greater 
Improvement  as  the  extent  of  the  convolution  Is  Increased.  For  Instance,  If  It 
were  desired  to  reproduce  the  function  over  192  halfwldths  using  four  decomposed 
functions,  a given  choice  would  give  a gain  of 


O'  1 + 12 


"T 

f 


192 


^«‘f6 


-t-  192 


6?’ 


192 

12 


16 


(43) 


It  must  be  emphasized  that  this  choice,  though  workab'^t,  la  not  unique  and  other 
choices  might  provide  greater  Improvement. 

The  algorltlim,  as  It  has  been  developed,  must  Include  protection  against 
lines  having  halfwldths  deviating  significantly  from  the  average.  For  anomalously 
wide  lines,  the  convolution  process  can  overflow  the  output  arrays  at  the  beginning 
and  ending  of  the  panels,  In  this  case  the  line  width  Is  reset  to  the  maximum  width 
that  can  be  treated;  twice  the  average  value.  A scries  of  + signs  and  the  line 
Identi  float  Ion  ave  outputted  when  such  a situation  U encountered.  In  general,  wide 
lines  are  not  a probleTi  since  the  widest  lines  do  not  vary  significantly  from  the 
mean,  This  esse  also  provides  protection  when  s sampll~  ' interval,  D'  is 
chosen  which  Is  not  oonslEtent  with  the  mean  halfwldth  oi  .he  lines  being  considered 
or  when  the  line  width  le  Incorrect  on  the  Input  data  file. 

There  are  a number  of  anomalously  narrow  lines  that  are  encountered  such  as 
the  high  J,  low  water  linee  of  type  measured  by  Eng  et  al.  For  such  lines, 
that  ks,  Oj  «:  a,  where  o^  Is  the  width  of  the  i'th  line  and  o Is  the  mean  width,  not 
only  does  the  algorithm  fail,  but  depending  on  the  phasing  of  the  sampling  Interval, 
tha  contribution  of  the  transition  can  be  missed  entirely.  Such  lines  have  been 

IT.  Eng,  R.S. , Kelley,  P.  L. , Mooradian,  A.,  Calawa,  A.B,,  and  Harmon,  T.C, 
(1973)  Tunable  laser  meaeurements  of  water  vapor  transitions  in  the 
vicinity  of  6 pm,  Chem,  Phys.  Letters  1B;524. 
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treated  by  setting  the  halfwldth  to  the  ttampllng  interval,  that  is,  - UV,  In  this 
case,  a series  of  — signs  and  the  line  identification  are  printed  out.  Although  the 
convolution  will  be  undereampled  for  such  lines,  information  pertaining  to  the  line 
la  retained.  The  only  strictly  correct  way  to  handle  these  narrow  lines,  is  to 
choose  a DV  equal  toa^/4  where  Is  characteristic  of  the  narrowest  line  encoun- 
tered. .Such  a choice  of  DV  may  necessitate  the  use  of  the  conventional  convolution 
method,  since  the  sampling  interval  U so  small.  Computational  requirements  may 
dictate  the  calculation  of  the  absorption  coefficient  over  a small  wavenumber  region 
around  the  lines  in  question.  In  general,  these  lines  are  relatively  weak,  and 
except  in  cases  such  as  laser  transmission  problems  the  method  used  hero  proves 
to  be  satisfactory, 


5.  EXTENSIONS  01'  THE  METHOD 

As  has  been  indicated  In  the  previous  section,  the  calculation  of  the  far  line 
wings  may  be  readily  incorporated  into  the  program.  However,  for  the  reasons 
stated  it  is  recommended  that  this  far  wing  contribution  be  handled  as  a continuum 
contribution  in  the  wavenumber  domain.  Similarly,  the  variation  of  line  shape 
with  molecular  species  may  be  incorporated  into  the  algorithm.  It  is  well  known 
that  the  wings  of  air  broadened  CO2  lines  are  "subLorentztan"  and  a very  slow 
function  specific  to  the  CO2  molecular  species  could  be  utilized  in  place  of  the 
Lorentz  very  slow  function.  However,  there  is  no  oompelllng  evidence  that  the 
subLorentzlan  behavior  occurs  within  6 oro*^  of  the  line  center  at  atmospheric 
pressure  so  again,  it  is  considered  more  appropriate  to  include  the  absorption  due 
to  the  far  wings  of  the  CO2  lines  ss  a properly  chosen  continuum  function.  The 
best  available  theories  indicate  that  the  subLorentzlan  behavior  is  wavenumber 
dependent.  Incorporation  of  this  effect  is  not  amenable  to  a direct  extenalon  of  the 
algorithm  alnoe  the  argument  of  the  Lorentz  function  is  in  terms  of  halfwldths. 

The  wing  modification  may  be  readily  added  to  the  calculation  by  Including  the  con- 
tribution Into  the  very  slow  array  before  the  three  arrays  are  combined  Into  the 
final  output  array.  Continuum  contributions  of  this  type  arc  Important  for  CO2  In 
apectral  regions  such  as  the  01001  band  center  where  the  wing  contribution  from 
the  Q branch  le  not  negligible  as  at  6tt7. 2 cm*^  and  beyond  the  band  head  from 
780  cm'*  to  800  om'^. 

The  convolution  method  deaorlbed  here  offers  even  more  significant  computa- 
tional savings  for  the  multilayer  case  as  It  le  generally  applied  to  atmoapherlc 
modeling.  The  sampling  Interval  la  determined  for  each  layer  based  on  the  average 
halfwldth  which  Is  a function  of  temperature  and  pressure  for  the  layer.  The 
remaining  problem  Is  to  merge  the  results  of  each  layer  using  appropriate 
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tnterpoUtlon  techniques  to  obtain  the  transmission  or  emission  for  the  case  being 
considered.  In  the  conventional  approach  to  the  problem,  tne  calculations  are 
performed  at  the  finest  sampling  Interval  required  and  is  the  same  for  all  layers; 
that  la,  the  sampling  Interval  required  to  perform  the  absorption  coefficient  calcu- 
lation at  60  km  la  the  same  ta  that  used  at  0 kin.  This  approach  requires  the 
calculation  of  many  more  output  points  than  is  necessary  at  altitudes  below  the 
original  altitude  with  a proportional  Increase  In  computational  effort.  With  the 
scheme  proposed  here,  the  computational  time  la  the  same  for  each  layer, 

b.  LINK  SHARKS  (ITilKR  THAN  LORKNTZIAN 

The  present  method  ia  clearly  not  restricted  to  any  particular  function  such 
as  the  Lorentz  function.  The  Doppler  line  shape  may  readily  be  tabulated  either 
directly  or  as  decomposed  funollons  depending  on  the  definition  of  halfwidth.  For 
many  atmospheric  modeling  cases,  the  line  shjipe  of  greatest  Interest  is  the 

1 O 

Voigt  shape.  Current  efforts  are  being  directed  toward  incorporating  the  Voigt 
line  shape  into  the  new  algorithm.  Preliminary  Indloatlona  are  that  the  Voigt 
profile,  computed  to  an  accuracy  similar  to  that  attained  for  the  Lorentz,  will 
require  only  slightly  more  computer  time  than  that  required  for  the  Lorentz  func- 
tion. Details  of  the  procedure  for  calculating  the  Voigt  shape  will  appear  in  a 
subsequent  publication.  Although  the  program  assumes  a symmetric  line  shape 
and  only  half  the  scanning  function  is  stored,  an  asymmetric  shape  nan  be  utilized 
by  storing  the  entire  function  and  making  minor  changes  to  the  program. 

7.  THKPHOOKAM 

The  program  has  t^en  written  to  uae  line  parameter  Input  data  consistent  with 
that  contained  on  the  AFQL  line  parameter  tape.  ^ The  line  data  ts  reformatted 
onto  a binary  file  which  contains  the  line  data  pertinent  to  the  molecules  and  wave- 
number  range  of  Intereat.  Thta  atep  has  been  taken  to  keep  read  time  consistent 
with  the  time  required  to  perform  the  oaloulattonal  part  of  the  program.  The  con- 
trol parameters  are  read  from  the  Input  file  and  written  to  the  output  file;  "tape  3" 
la  the  binary  file  containing  the  line  parameter  data;  and  KKFILE  la  the  binary 
output  file.  KKFILB  contains  a header  record  which  Includes  the  Identification 
Information,  SECANT,  temperature,  pressure,  molecular  Identification,  and 
molecular  column  densities  of  the  homogeneous  layer.  The  first  record  for  each 
output  panel  Is  a header  record  for  the  panel  which  contains  the  wavenumber  valuea 

18.  Voigt,  S.  MOnoh,  Ber.  (191!))  p.  d03. 


of  the  firat  and  last  absorption  coefficient  values  of  the  panel,  the  wavenumber 
increment  between  output  points  and  the  number  of  output  points.  The  second 
record  of  the  panel  contains  the  array  of  absorption  coefficient  values  resulting 
from  the  convolution  oalculation.  The  current  version  of  the  program  outputs  a 
maximum  of  3400  values;  in  general,  the  firat  and  last  panel  are  shorter. 

The  program  consists  of  the  main  program  HIRACC;  subroutines  SHAPE, 
MOLEC,  RDPILE,  CONVFN,  and  PANEL;  and  the  function  QVRFAC.  The  overall 
strategy  of  the  program  la  Indicated  in  Figure  0.  All  the  subroutines  are  called 
from  the  main  program  and  the  flow  of  the  program  la  easily  traced.  Subroutine 
SHAPE  sets  up  the  three  convolution  functiona  XF,  X8,  XV.S  used  to  define  the 
Lorentz  function  from  0 to  48  helfwtdths. 

Subroutine  MOLEC  in  conjunction  with  function  QVRFAC,  makes  the  molecular 
identlfloations  with  tha  line  parameter  file,  and  determines  the  correction  factors 
for  the  line  intenaltlea  <SCOR)  and  the  halfwldtha  (ALFCOR).  Tha  quantity,  SCOR, 
is  the  correction  factor  due  to  the  temperature  dependence  of  the  vibrational  and 
rotational  partition  sumn.  The  vibrational  partition  aum  is  calculated  for  a given 
molecular  type  aa 


Qv(T)  - n 


X 


t>l 


«4) 


where  Is  a fundamental  vibrational  frequency  and  d^  la  the  degeneracy  of  the 
vibration.  The  temperature  dependence  of  the  rotational  partition  aum  is  given  by 


<40) 


where  F ■ 1 for  linear  molecules  and  1.  S for  nonlinear  moleculea.  The  reference 
temperature,  T^.  la  taken  aa  390  K,  conaUtent  with  tha  AFQL  Line  Lilting.  For 
further  dlicuaston  of  Ihaae  topice,  eee  Herzberg^^  pp  S03  It.  The  partition  aum 
caloulatlona  are  performed  by  QVRFAC  and  the  neoeaaary  moleoulur  parameters 
are  contained  in  data  statements  In  subroutine  MOLEC,  The  quantity,  A LFCOR, 
la  tha  oorraotlon  factor  due  to  the  praasura  and  tamparatura  dependence  of  the 
oolllaton  broadened  halfwldth.  The  temperature  dependence  of  the  halfwldth  has 
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Figure  6.  Flow  Diagram  for  HIRACC  Program 
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n 

been  taken  an  (T^/T)  ' although  calculatlona  baaed  on  the  Anderaon  - Taao  - 

Curnutte  theory  are  reported  to  give  a temperature  dependence  « (T  /T)^‘ 

20  ° 

(Varanaal^”). 

Subroutine  RDFILE  reads  the  blocked  binary  line  parameters  over  the  wave- 
number range  for  which  line  data  la  required.  The  lino  parameters  Include  the 
wavenumber  value  of  the  tranaltton  {QNU,  cm'M.  the  Intensity  of  the  transition  at 
208^K  (S,  cm' V(>nol/cm^)),  the  collision  broadensd  halfwtdth  at  half  maximum 
fur  886^K  aind  atm  pressure  (ALFAO,  cm*^),  the  lower  state  energy  of  the  transi- 
tion <EPP,  cm'M,  and  the  molecule  Identification  number  (MOL).  If  the  line 
parameter  data  Is  Insufficient  to  complete  the  specified  calculation,  the  message 
"end  of  file  on  disk"  Is  printed  on  the  output  file.  If  no  further  line  data  Is  required 
IDATA  Is  set  to  1,  and  control  Is  returned  to  the  main  program. 

At  this  stage  of  the  main  program,  an  effective  optical  depth  Is  calculated  for 
each  line  which  la  dependent  on  the  column  density  of  the  layer,  the  secant  of  the 
angle  through  the  layer,  the  temperature  of  the  layer,  and  the  halfwtdth  of  the 
line,  olALFI), 

The  effective  depth,  u', 

(46) 

where  w Is  the  absorber  column  density,  E"  Is  the  lower  state  energy, 

Qyp  * quantities  have  been  previously  defined.  In  terms  of 

the  program  coding  the  affective  depth  appears  asi 

aFDPTH*RECALF'>W*SEC*SCOR*EXP<EPP/XKTFAC)*(l.  -BXP<-ONU/TEMPO))/ 
(1.  -EXPl-aNU/TAVE)). 

As  prsvlously  discussed,  the  proper  sampling  Intsrval,  DV,  should  bs  0. 38  times 
ths  average  line  halfwldth.  If  the  halfwldth,  ALFl,  la  leas  than  ths  aampltng  Inter- 
val, the  halfwldth  la  aat  to  the  sampling  Intsrval  and  a seriaa  of  mlnua  algna  la 
wrlttan  to  ths  output  file.  If  the  halfwtdth  excaeda  a maximum  value  (ALFMAX) 
where  ALFMAX  • BOUNO/40  and  BOUND  la  the  maximum  value  In  wavenumbers 
over  which  a line  can  be  calculated,  the  halfwldth  la  resat  to  ALFMAX,  and  a sarlaa 
of  + signs  Is  written  to  the  output  file,  bioluded  In  the  records  Indloatliig  the 


30.  Varanasi,  P. , end  Ko,  F.  K.  (1077)  Intensity  and  transmission  measurements 
In  ths  Vfl  fundamental  of  NeO  at  low  temperatures.  Thirty  Second  Svmnaalur 
on  Molecular  Spectroscopy,  Paper  RF»5.  Columbiial  taO. 
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r«a«ttlng  ot  the  halfwldth  la  the  wavenumber  value  of  the  traneltion  (QNU),  the 
Intenatty  (S)  and  the  halfwldth  (ALKAO)  values  of  the  transition  from  the  line 
parameter  tape,  the  calculated  value  of  the  halfwidth  (AI.1FI),  the  value  to  which 
the  halfwldth  (DV  or  AIjKMAX),  and  molecular  Identification  number  (M),  If  the 
number  ot  halfwldth  changes  (NCHNQ)  exceed!  100,  the  computation  le  terminated. 

Subroutine  CONVFN  la  a tightly  written  subprogram  in  which  considerable 
effort  has  been  taken  to  mlnlmlxo  operations  in  the  DO  30  loop.  This  subroutine 
performs  the  triple  convolution  of  XF,  XS,  and  XV vS  with  the  line  data  putting  the 
results  in  the  proper  elements  of  FF,  SF,  and  VSF  respectively.  A simplified 
flow  diagram  appears  in  Figure  7.  Control  indicator  IPANEL  is  set  to  IDATA  If 
the  DO  loop  over  the  lines  (40)  la  selisfled  Indloetlng  whether  a panel  is  complete 
or  more  Unas  are  required.  If  the  line  DO  loop  (40)  la  not  completed,  IPANEL  Is 
set  to  1 Indicating  that  e complete  penal  has  betn  calculated.  Control  la  returned 
to  the  main  program. 

If  IPANEL  has  been  set  to  1,  aubroutlne  PANEL  is  called.  Subroutine  PANEL 
performs  a four  point  Lagrange  interpolation  of  the  VSF  array  into  the  SF  array 
and  the  SF  array  into  the  FF  array,  thua  oomblnlng  the  reeults  of  three  Independent 
convolutlona  Into  e final  result.  A gtnsral  flow  diagram  of  panel  is  given  In  Fig- 
ure B.  Care  is  taken  to  etora  array  values  required  for  the  Interpolation  of  aub- 
eequent  panels.  VFT  Is  the  wavenumber  value  of  the  first  element  of  the  FF  array, 
which  is  common  to  the  first  elemmt  of  tha  SF  and  VSF  arrays.  A binary  header 
record  is  written  to  the  binary  file  (KKFILE)  for  each  panel  which  Includes  the 
wavenumber  value  (VIP)  of  the  first  element  of  the  panel  (FF(NLO)),  the  wavenum- 
ber value  (V2P)  of  the  last  element  of  the  panel  (FF(NHI)),  the  wavenumber  Incre- 
ment (DV),  and  tha  number  of  absorption  coefficient  valuta  outputted  (NLIM).  Tha 
second  binary'  record  contains  the  NLIM  valuea  of  the  absorption  coefflotsnt  from 
the  FF  array.  The  arrays  are  appropriately  shifted  and  reset  In  preparation  for 
tha  computation  of  lubsequent  ponela.  Control  is  again  relumed  to  the  main  pro- 
gram, HIRACC. 

After  a panel  Is  written  to  KKFILE,  e record  is  written  to  the  output  file  indi- 
cating tha  current  value  of  the  time,  the  time  spent  In  ROFILE,  in  CONVFN.  and 
In  PANEL  (the  units  are  aeconda).  Also  Included  In  this  record  ere  tha  firat  and 
lost  wavenumber  voluaa  of  the  panel,  the  wevenumtier  Increment,  and  the  number 
of  values  In  the  panel.  A aeoond  record  is  written  to  the  output  file  indicating  the 
average  value  of  tha  holNrldth,  the  number  of  llnea  reed  since  tha  lest  panel  waa 
completed,  and  the  total  number  ot  lines  read  since  the  Initiation  ot  the  convolution 
oaloulatioB.  Control  la  returned  to  statement  10  If  the  calculation  la  complete,  or 
to  statement  140  to  continue  the  convolution  prooeee. 
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Figure  7.  Flow  DUgrun  for  SUBROirrmiC  CONVFN 
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Appendix  A 


A llutlng  of  th«  Fortran  program  HIRACC  la  gtvan  In  TabU  Al,  togethar  with 
tha  loadar  map,  the  Input  data,  and  the  program  output  for  the  reaulta  which  appear 
In  Figure  6.  in  the  output.  It  ahould  be  noted  that  one  Una  la  encountered  that  haa 
a halfwldth  which  la  larger  than  that  allowed  by  the  program,  Thla  line  at 
3fiSB.4B30  cm'^  la  an  HDO  line,  molecular  type  1,  which  haa  an  erroneoue  half' 
width  on  the  AFQL  line  tape  uacd.  Tha  width  for  thla  Una  la  raaet  to 
ALFMAX  ■ 0. 08  cm’^.  Two  HjO  llnea  occur  at  3979. 0810  om*^  which  have  ano- 
maloualy  narrow  halfwldtha.  The  valuea  of  theae  halfwldtha  are  correct  and  are 
the  high  J,  low  water  llnea  of  the  type  meaaured  by  Eng  at  al,  ^ Tha  wldtha  for 
theae  llnea  are  react  to  DV  ■ 0,01  cm'^.  The  beginning  and  ending  aeotlona  of  the 
output  panela  are  printed  out  Including  the  location  In  the  FF  array,  the  dlmanalon- 
leas  abaorptlon  coefficient,  and  the  wavenumber  value  aaaoctatad  with  the  respec- 
tive elamenta  of  the  FF  array. 


1.  Eng,  R,  S. , Kelley,  P,L.,  Mooradlan,  A.,  Calawa,  A.R. , and  Harmon,  T,C. 
(1073)  Tunable  laser  meeeurementa  of  water  vapor  tranaltlona  In  the  vicinity 
of  8 itm,  Chem.  Phya.  Lettera  10i534. 
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Table  A 1 


LlBtlng  of  HIRACC  Program  (Cent. ) 


ii* 

Ut 

ire 

ITS 

Its 

111 

111 

lit- 

til 

IIS 

III 


MIRICC 


til 


ill 


111 


111 

0 


e 

111 

c 


V 

Ul 

c 

c 

c 

III 

«» 

SIS 

SIS 

111 

SIS 

StI 

SIS 

111 

SIS 

SSI 

SSE 

SIS 

SSI 

sss 

SSI 


IT  ItrOSTNIIt.ll.l.t  CO  TO  ill 
lONTalCNTTl 

liTIalLTIKII'lLfOOKM) 

IT  IILTI.M.DSl  00  TO  111 

TRINT  SSS,  SIM<n,SiII,*Lr*OIIttli.rt,DV|M 

ILTIaOV 

NCNNCaNCHNCTI 

IT  IILTl.LT.  ILTkail)  CO  TO  ItO 

MINT  SSI,  6NUcn,S(ll,*LTA0(II,*LTl,ll.rHIII,H 

ILTIalLTHTI 

NCMNCaNeHNCtl 

CONTINUC 

SUNU.T>IUHiLT«HTI 

noiiTU)>u/iLTi 

tT0TtM(Il«|F0TrH(I>»£«N«i:M«It*l(t<tT*0|T«C»LTCI>> 

I 1 1.  -t  *H  -ONOt  I)  /T»«t ) / a ,-**T(-6NUII  > TTINNO)  I 
nONTINUt 

IT  (NnNN6.CT.lQII  CO  TO  l«3 
CONTINUr 

GILL  OONSTN  <GHUtrTO»rH,S(CALT,NOi,TT,IT,V|T,l(T,l(S,«VS» 


IT  (ITANtL.eO.Q)  CO  TO  101 
OONTINUT 

OALL  PANEL  (TT,CT,«Sr,i(KriLE) 

CALL  irCONO  IfINEI 

PNINT  SSS,  TINC,TIM»0T,TINCNI,TINPHL,V1P,V1P,0S,NLIM 

ASALPTlUNILF/FLOATItCNTl 

LINCNKLINCNTtICNT 

PRINT  Sll,  AVALTiICHT.LXNCMT 

SUNALTtI. 

ICNT-I 

IT  (tSTOP.NC.tl  CO  TO  LSI 
CNO  TILC  NKTILC 
CO  TO  II 

CONTINUE 

tfiJP 


rOANAT  (IHtI 
fOCNATIIl.riil) 

TOIMArilN.TAtll 

rORNATlIS.TtS.CI 

P0RNAT(»0  OUTPUT  TILS  NO,  ••  III 
rORNATI*!  SECANT  TllaSI 

fORNAT  (tl«IHll,llt  .. 

rORNATITI  PRCSKNR)  ••  Tit, ST  •!  TENPdO  •*  Pll.tl 
rORHATdSNI  OOLUNN  OCNIirriNOLECULHTONT*?)  , 

1 T /(SK.AI,  * ■ * IPIll.l)  I 

TORNAT(»|^oJ(0N-tl  • •fll,N/T|  VllCH-l)  • TTltiST*!  VtCCN-ll  ■ • 
*fORN*TUTHO  I'lOUNOlCM-ll  • ,TI.I.I 

f diMAT  •Alltel  I • If  It  1*9* 

FOVNAT  If  It 

TORNAt  (INI,  llK,*TIHtT,llN,*REA0»,S»,»WMy0lUT16NP,llN,*PANCl*,SN, 
l«Sl*,llll,*Vt*,llT|*W*CN|*NSALUES*,T()Tll,l,ITiS,S,Tll,  A,  III) 

TORNAT  (*1  AVI  RICA  NIOTM  • •Tll.S,*  NO.  LIHIt  ■ *Illi 
t » TOTAl  NO,  OT  LINIC"»III/T» 

ENO 


It  till 
lOlllO 
lOtilO 
lOllfO 
lOlINC 
OCIASO 

lullsi 

CtllTI 
00  HID 
OOilSO 
001710 
00171Q 
OCIIIO 
001710 
061710 
001710 
001760 
061770 
0017AO 
0017SO 

ooliio 
001110 
OOlltO 
0011 10 
ObliSO 

00  use 
OOIISI 

001170 

ooino 
00 1100 
00  ISIS 
OOISIO 
OOlSIO 
061SIO 
lOlSSO 
OOISIO 
lOISOU 
00IS70 
OOISIO 
OOlSSO 
OOIIOO 

ooini 

OOtlZO 

OOtllO 

OOtlAO 

OOIISU 

OOtliO 

OOtITO 

301100 

OltlSO 

OOlllO 

OOtllO 

ootito 

OOtllO 

RHINO 

IIISO 

OOIICO 

OOtllO 

OOtlM 

ootino 

oottno 

autiio 

oottto 

OOtllO 

OltINO 


til 


Table  Al.  Listing  of  HIHACC  Program  (Cont. ) 

. • F’N  11.5 

suBnoutmi  sMift  *>*'’•* 


urni-ni 


1 


It 


tt 


It 


ft 


31 


It 


tl 


tt 


to 


It 


It 


M 


tl 


II 


II 


vutwut  iNt  tuiPf  (»ri  »s,)<ist 
'nitfMSioN  yriii,  «s(t>i  tvtdi 

KtCtNf  ( Hi)  •i*/«  i**«f 

)!i(Si”*tt*i’*j 

on  UMl* /<!•♦*•••*•*** 

HAT  CM  AT  1 1 HALfHIOTMI 
I l*iii 
Al>ll<Ill 
Hl'IKIl) 

Ql«Ql  (tl> 

1 HATOH  AT  It  HALFHIOTHf 
AfoAKItl 

at>io«i» 

CfaClOt) 

TOTAlali 

00  11  ]<t«NfNAI 

KKU  •IICHl'l  flO****  ••  ••  • 

suN>ir(ti 

no  tl  JJ'fiNF 
T•rL0IT(JJ-H•01<F 

lUNHtUHtXf <3J> 'll 
tl  CONTINUf 
KFINF)«I> 
lUHaflUHtOKF 
TOIAt«TOTAHlUH 

00  31  1«1«NSHAX 

**  5»tt>'Fio»l»<*F*<l*****^''***** 

UllHfJoiT  (HHrWai(S*l.lll 

01  Al  J^liNSl  . 
fafLOIT  UJ*H*D»* 

in JjNwCK*  him  <« I -VIFM 1 «» » 

IM.IUHtXKJJI  •!. 

Al  OONTIWI 

mtlaMITl 
00  II  JJaNlte>HI 
(•IVUIT  >*<■«'> 

xil  JjJallOM*  IXlOIMICXtl'Wm***  • 

IUH«IUHaKI(Jtll*t> 

II  CONTIHUf 

KKHtlAla 
lUNatUAtDIS 

totaktothtiw 
00  II  taliNIIMIK 

*'  mliilSioMAWewt* 

lUNalliai 

nvii«ho*t  (MM*i/axv««i.iia 
00  tl  Ji»liN»ll 
XaFI.0AT4JJ«li*0XV« 

yf  ■ If  M 

tUN*IUMTKftUJ>*l< 
tl  OONTIHUC 

NHIAiHllllI 

OO  11 

««ri.oiT<JJ<i>*Oxil 

IlKJJlalieiXTIXtOIMK*** 
igMafUNTXIII 
It  eONTlNUI 

IUH«IU*<*(K»» 

TOVIAOtOTIlTIVI 

MTum 


oottto 

Jt.23«0 

1022/1) 

0022<0 

0C22tO 

01.2.100 

002110 

002120 

002110 

0023it0 

0021I0 

lOtllO 

0021/I 

OOtlll) 

DOtltO 

002AOO 

lOtAlO 

OOtAIO 

OOtASO 

OOIAAO 

OOtAIO 

OOIAIO 

OOtA/O 

OOO'AaO 

OOtAIO 

oottio 

uottto 

002120 

002110 

OOtlAI 

OCItIO 

oil  till 

0021  to 

00lt|l 

oottte 

OOlAOd 

lOlllO 

002ltl 

OOtllD 

oAtlAO 

OAtni 

eitiii 

dOiiri 

OOtllO 

OOIAIO 

ootrio 

oitrii 

iiirii 

OOIUI 

I0I2AI 

lOITII 

OOttAI 

lOtT/l 

lOttlO 

DOtrii 

OOtlll 

lOtllO 


lOtlZI 

001130 

OttlAO 

OOtllA 

OOtlll 

lltlio 

OOtlll 

OOtlll 

OOtllO 

lillio 

lOfItO 

OOtlll 

OOtlAI 

OOtllO 

OOtll) 

oiiiro 

oonio 

OOtlll 

OOlllO 

tlllll 

iiiuio 

•AlUIO 

IIIIAO 

ictaio 


4i 


Table  A 1.  Llatlng  of  HIRACC  Program  (Cent. ) 


sun  ROUT  I N(  HOLCO 


RTN  It.Rtlilli 


8F/11/TT 


SUIROUtlNf  HOLSe  (INOimeLlOiMtKIllO.TCHM*  tIMR»RO|R(5COII|»lRCORt 
OIMNttON  HNOLIDUl  iSOOHdt  «RLFCOR(li 

DIRINSIOM  M<r|til|NO«r«lit  

OOHHON  /RMOtM/  HH  « Jl  tM»«  FI  . Rl«  fl  iNi  < Tt  t H*  ( f ) I F I i 
1 HI  (Tl  <N3iri  iNltCFI  •QVtiri  iROTriCIFI 

OITR  MHUtiNRtll  ,l(HlliNHlT»Rtt  lltWfU)  |RSI  i(  iMSH)  )mi<  1)  tMHIl)  t 

t OVdH.ROTrRCIH/  , . , , , , 

I IH  Mil  ill  SSIliTiti  ISSSiSili  IFSSilili  0i  iDt  1 ■dp'** . * ’ 

tllTI  HM(FI  iSVOt  iMl  (It  iNllt)  iNtIt)  |Mf  (I)  |HS  t*  > itJ  <*  t iRRIIl  |Wk  tl  I l 

t niQIIliROTHOdl/  „ . 

I tH  cot  ill  IJRkiOili  »(tF.li*i  IlkSiIili  0«  iti  liQlkSi  I • 0 / 

OITA  HRISlikVdl  imtlliMHItiRf  «IlittHJI|Rt<ItikH*t[UktJliNktJli 
j qvo(iiiROTr*cni/  ^ ...  .....  . . . 

X IM  Ol  liolilill  TtQ.klll  lOklilili  0i  lOi  lilkkki  Lik  F 

OITA  MR  tkt  iMV»kt  lUlIkliMt  <kl  iltIlkliNI(kliRJ(k  liMHk)  |UkUI  |Nk(kl  , 

A OVO(kl iROTFAOTkt/ 

k AH  MO  ill  KAk.tili  klAiAiIi  IIIliAili  Oi  |0i  l.llAFi  liO  / 

data  HHUliHSdl  |Hl(ll|MHSI|H»IAtiNIITtiHJ(l)lHJt»t|Hk(»l  iNklk)  i 

A a»o(S)iROTrAO(si/  , , 

« AM  CO  ill  Ilkliliti  Ml  iHi  l.AtIi)i  LiU  / 

Utf  A MH«*liHV<AI|Hl(AliHt(»l  iHHAt  |M*  (Al|HJ(*)iHIUt  iHklll  |Nk(AI  i 

A AM  CMk  iki  ISlki  lii  liSIilifi  SOAiiSiSi  lIAFiSiIi  l.AOAAi  Sif  / 
*OA1A  HM«r!iMV<Fl  |Hl(FI|Nt«TliW(»liMHFl|MH  TtiNUFI  iMktFt  iMktri  i 
f OVOdliMTFAOATI/  ...  ....  , . j 

f AM  01  ill  tkl.TiAili  «•  ill  ••  1*1  ••  I*'  iitl*Fi  liO  t 

MOklC  MAMI8  TMC  NOKOULAR  tOINTirtORTlONt 

ALROOR  IS  TMI  FACTOR  « RMICM  TMt  COUIltON  MIOTH  MUST  Rl  OHANBID 
out  TO  THS  oeClMDtNCI  ON 

TMI  TtMRtNAIUNt  OeKKBtNOI  IS  TAKIN  AS  (TO/TI»»S.S 

SOON  IS  TMC  factor  RT  tRtlOM  IMS  ItNt  INTINSITT  IS  OMANMO  OUI  To 
TIMFIRATORf  OtWHOfNOI  OF  TMt  VI*  AMO  ROT  FARIITION  SUMS 

IF  IINOiCS.il  SO  TO  fl 

NOIHhT 

NVOIMik 

OU  IS  N«liNNOl(0 

N(Nill>Ht|MI 

HIN|>I«MI«NI 

MINill«HS(NI 

NINlkl«Nk(MI 

ND<Nitl«NI«NI 

NDINill>NIINI 

NOlHill'NMRI 

NDIN,kl>NklNF 

HNOtlOINlHMNIMI 

RCTWN 

CONTINUI 

RRATInaFFRO 

TRATlO»TtMFS/TrMF 

SeOR*|N»R8iRFAO(M|TtNF|TRAtIOiOyl|ROTFA6iM.NDiMOIMiMVO|MI 

AlFeORIN|iFRATIO»«TRATIO»*S.»t 

CONTINUr 

RFC  URN 

CNO 


scioso 

SAIAFO 

OOlSRO 

00J04I 

OIILBO 

003110 

001120 

D01130 

eoiikc 

OullRO 

0011*0 

OOIIFO 

OOIIH 

001140 

OOllilO 

OOlllO 

00122(1 

0112JO 

OOlfkO 

0012*0 

0012*0 

OOI2FO 

001200 

00) 240 
OOISSO 

01) 110 
011)20 
0011)0 
lIDkO 
00)110 
•l))*l 
OOMFI 
OOISIO 
001)40 
OOlkUO 
DOlktO 
OOlkOO 
OOlklO 
OllkkO 
lOlkSO 
Ollk*(l 
IDkFI 
OOlklO 

ooikse 

01)110 

001*10 

OOlSlO 

il(|]*)l 

OOltkO 

01)411 

001**0 

oii*ro 

00)410 

011**0 

003*00 

001*10 

001*20 

011*10 

0Dl*k) 

Oa)**l 


r 


TkbleAl.  Listing  of  HIRACC  ProgrMn  (Cont.) 

' I 


ruNcrtoN  ivitrte  Fii/i’k  oft»i 


rfN  4it*kllt  •7/11/7/ 


1 


II 


• UNCTION  0»»rAOIM,TIN7|TNIT10tO¥liWT7*0«N|NO»HUIN|NVOXN» 

OtHFHKlON  OVOlllfMTriCIlt 

OINfNItON  N(NO|H«NVDtMI  ,HO (NDIN,N«DtHI 

oval. 

DO  11  laliNVOtN 

17  INOIN.  11.(0.01  SO  TO  (0 

tVal.-tX^I-NlH.ll/TIMNI 

!•  (N0III.0T.il  tVtlV*«ND(H.II 

to  ov«av/tv 

tl  0VN7*0a  lOVO  (HI  /OV  (•  I ( TMT  I01»»  WTMC  (HI  I 

OITUNN 

0 

(NO 


IIISOO 

tijsri 

00J4OI 

OSIMI 

iijroi 

•13711 

• I17II 
tgirii 

• DI7U 
•I17H 

• el7tl 
III77I 

• tl7IM 


tUSNOUTINC  NOriil 


rt./7l>  ONTal 


FTN  •.•••ll. 


•r/ii/7T 


II 


tl 


31 


M 


•I 


•t 


70 


m.lllTI.IITOO  . ____  


til 


1(11.  AkrlKlI.  CNNdl.  NOlltl 


OININIION  INUltl) 
lOITKI 

UlO  (11  VHtN.VMIX.NNCe 
IF  ItOFlin  ll.ll 
IF  (VMK.OI.SIOTI  00  TO  SI 
HIO  (II  VMtN 

NIIO*'lil  (INUdl  .1(11  (ItFIK  1,1. ((►1*111  .NOLIII.I»l.NN(0> 

UOat 

INI*N«0  , _ 

IF  IVMIN.il. VOOTI  00  TO  II 

00  II  lal.NNIO 

ILONi 

IF  lINUIII.SliVMTI  10  to  II 
OONTIMUV 

IF  (VNIX.LC.VTONI  NITWN 

00  II  lai.ime 

IF  (INUID.IT.VTOFI  10  TO  7| 

OOHTINUC 

IF  (IMt.lf.HNIOi  leiTI't 

ItTIMIN 

FNINT  Ml 

II0F«1 

KTUNN 


«l  FOIHIT  (•  (NO  OF  FlU  OH 


•t.irti 

••soil 

iijiii 

IIMII 
11*111 
IISIII 
lit  — 
III- 
11*11 
•1*11 
•1*11 
II 101 
II**' 
Ml 
Ml 
Ml 
MMI 
MMI 
II  III 
II  »< 
MM* 
MMI 
•IMI 

•I  Ml 

•IMI 
MMI 
MM< 
IIMI 
IIMI 
IINW 
IIM' 
IIMI 


Table  A l 


L-Utlng  of  HIRACC  Program  (Cont. ) 


SUIHaunHC  CDNVrN  OI»t«l 


TTN  It.'in.llt 


i 


« 


It 


to 


It 


It 


*1 


«tt 


(I 


ft 


•I 


ft 


I'l 


n 


It 


c 

c 

0 

c 

c 

0 

c 

c 

n 

c 

0 

n 

0 

C 

e 

e 

c 

0 

Q 

C 


It 


>1 

kl 

e 


It 

ft 


fuifouTiMF  ooNvrN  ieNUtFFtFfM,aio*cr,HOi.i)'F)tF<viriiF,ratitvi> 
OOMFUtl  /tsiwy  aiOF>iltHINilLO.tHt.VIDT>»TCIFitliftiOVtl«rTi  tCOr,IF«N 
ItltlOttlilltTO* 

OOHHQN  /Mtx/  NMriaKritlF|)WS(OI(fiNt|NNVI)OIIVttNI(Vl,NrN»ll,NIHAX|NfSN 
lAK 

OOHMON  ^tuai/  MAlir>M*«f  •H«l(VS»NLttr>NLlNltNltNUStHLOiNHI  t0vs«0tvs 
COHMOK  /XlIHt/  TUffliTtHROFtTtMlHVirtHFNL 
OIHtMlION  GNUUIt  CFOFtHIll.  *tC*LF(l»i  FFMI*  SFIll 
DIMIMtlON  NOLU) 

OIMINllOH  Xfdii  XlUlt  xniltt  VSFIll 


TM  f(Fy  tlOH  FUlCTtON  IS  UttO  TO  DtIVC  TMI  ODNVnLUTtON  PKOCISS 

IHI  «t«V  II.4H  CONVOLUTION  IS  PCNrOtNCD  TNOH  •NMVS  TO  «)<MVS  IN 
HILTHICTHS 

OXVt  It  THI  INKNVAL  » NNIOH  TNC  VINT  ILOM  FUNCTION  IS  ITONtD  IN 
HtLFNirTNS 

JJ  IS  TNI  LlOtTION  IN  THC  VtF  tPNAT  FOR  MNICH  TNI  CDNTNItUTION 
TO  INC  ANSONFTION  OOCFFIOItNT  FROM  THI  PTH  LtNI  IS  CtLOULtllO 

Ills  It  TMF  tRRROFRUTI  INOCN  OF  THI  XVI  FUNCTION 

VFT  It  INI  UtVFNUNRU  VALUF  OF  LOOM  ION  1 OF  THC  FF,SF  AND  VST 
IRFITS 

OCFTMIa  CFFIOTIVI  tnORRTION  CONFFIOHNT  OF  THf  PTN  LINI 


euL  MooHo  piNini 
RtTVX«OVVI/OXV1 
FVSaFLOtT  INNVf  MOXVt 
flONF>QVVt/aV>t« 

OONS«OVVl/OVt*t> 

IF  <IL0;CT,|HII  to  to  ft 
on  Nl  MILOiIMI 
OlFTHIalFOFTNUI 
IF  (DIFINlAC.Itl  GO  TO  Nl 
tSiORORCCILFI  tPRITVX 
tIRT«IINU<lt-VFIt/OVVS 
IFV>>F«VI&ORC 
JN*tN4TIHT*tFVt>»t. 

IF  <<IMI.LC.H*VXSI  60  TO  19 

ILIST*I-i 

60  TO  It 

JMIN>«nNT>NFVSl  *F. 

JNINFallXNT'OOHri 

9HINSa1IINI«OONI> 

IVIalFlOIKJNIN-lt’IlNTPISLOFI 

ll■FVS•(FLO•T<  JNINSl'IINTaoORSPZSLOFI 

ZFilVIT  IFLOIT  UNINFI>ZINTaCONF|»ZtlOFr 

OU  Jl  JJaJNINf JHIX 

JSaJNINSFJJ 

JF>JN|NF*JJ 

ZVfaZVKZILOFI 

ZtaZSilSLOn 

IFaZFTlILOFf 

nVRaMSIIVIIFI.f 

IZ|Htlt(ZII*L,t 

iiF>*itirFi*t.i 

VSF<9JI>«IFI9JI  tOK*THI*XVt(IIVI) 
tF(JS|iSF(JS|*BIRTNI»taUZ« 

FF<  jn>FFI9FI*0IRTNt«XPIIIFI 

CONTINUI 

OOHTINUI 

HIST  MINI 

ISOIIae  IF  ROK  Olfl  RI4UINIO  ••  lOlltat  IF  NO  NODI  OATI  RtOUINC 
IRAWlalBITC 
•0  TO  II 
eONfINUt 

IFIWUI 

iiilL^tttoRl  <TtN«l 
T IROHVaTIHCNVNT  IRC'TINM 
RITURN 


(NO 


OONIiO 
OONIZO 
OINIIO 
OUNUO 
OONitO 
004110 
0041  FI 
104110 
0041t| 
004100 
004110 
OONIZO 
004110 
004140 
UONlfO 
004160 
004110 
004II0 
044100 
004100 
I04I1I 
004)10 
904110 
01.'  4140 
004100 
004160 
oi4iro 
004)11 
104)00 
004400 
0I44II 
004410 
0044)0 
II444I 
014400 
004460 
004470 
014411 
014400 
014001 
014010 
004ICO 
004010 
904I4U 
004100 
114061 
904IFI 
004000 
004000 
004600 
004610 
9046ZO 
0046)0 
904640 
90 4600 
004660 
004670 
004601 
004600 
904700 
004710 
004710 
004710 
004760 
II470I 
014710 
014770 
004770 
104710 
004709 
004011 
004119 
114000 
9040)0 
014161 
014000 
014060 
0049  70 
II4IOO 
004000 
004000 


0 

0 

0 


44 


V'.‘ 


Tabla  Al.  Lilting  of  HIRACC  Program  (Cont. ) 

SUtHOlUlNi  PtMll  Tit/Tl.  OPTel  rTN  il.S  + lim 


SURNOUTINt  PANEL  CPtAriVNrtKKrtLEt 

COMMON  /XSU1/  NFCPIiLIHtNitLOt1HI|VR0r,VTaP,VliV7iI]V,vrl,  tf UP  (IMAM 
I'Ll  IDIIAiMToP 

COMMON  yRUtl/  HANr|HAIIS|NAI(V^,NLINr.NLtH$,NlIMV^,NLO|HHI  iDVSiDWVA 

COMMON  /IMAMCL/  « AM| VtPiML IMiNiHt'T « NPT S 

COMMON  i>«TtNt/  TIME  ,nHllorittNCN«iTllMMM. 

niNINIION  PMItl,  iruii  Nsrill 

DIMENltON  NMITUMdl 

COLL  KCONO  ITtMftI 

NIlMt. /If*. 

«IIn|i/1Ai 

J8T8M«C 

l|NHl'(«t-VPn/OV*l.f 

IP  INNl.et.NHMtl  IllOP'l 

IP  ll«TOPi(Qil)  NMIaNNHI 

JMIIPMNI.O<NaTI 

JjMiNMI'NPTt 

LINia«<NLO>ll/N-l 

LIHMiNlAMIAttt 

ILOH  PUNCnON  EALUCt  NCIOIO  PON  SUMIOUINT  PAHELt  ANP  SAELD 
DO  tl  J>l«i 

|PITON<JMtPUIMMI*J-tl 

on  n j«iiiNi.OiiiNMi,A 

JEI>(J<1I/A*1 

IPCJI>«Pt4l'EfP<Jtlt 

lPIJA|l>IP«JMI«AE«<»tPOES>tl»E01*ElPIJE«>*EOC*EJPIJEI»l)*>AI*ElP 

IdEEtEI 

tPI  J«tl  ■ir(J'fM«tt*«EIP«JEI«ll«ElP|JEi*>l»»EII*  lElP  IJEIMEIP  IJEt* 
Ull 

fPUAlt  'SPI  J*]»  ««El«VtP(  JE«-tl«>OP>ESP  IJESIaKII'VAPIJES' tl'XOO'ESP 
XJElAEI 
OONttNUI 

no  10  J'NLOiNMIiA 
JU  U-tiPA*l 
PP(JI*PP<JI*IP|Jt) 

PP«J*ll«PP(J*ll*«OI«IP«OS*l|AHI*lPUtl  AAIII*tPIJS'll  AEOJ'tPUS'II 
ppiOACi  APPiJtii  *«uA«tr<Jt>i»*fP  us»r)i»EtiE(«PiJiiASPU>ni) 
PPlJAltAPPU'll  •EII*IPU>>ll*A0tPlP|J*l<X01MP(,tl»ll  'EIS*IP(  Jl'l) 
COMIINUI 

IP  IN»t«,IO.OI  eo  TO  AO 
00  A|  jANlOiJNEPN 
NNnVfrlAPLOATOItEOV 
NNINT  tilt  JiPPdhVP 
00  00  JNJJPNillMI 
E|MEPT*PlOAE(J«tl’aV 
POINT  fill  JtPPUliVP 
eONTINUI 
00  TO  j>l)0 

tPIltWiAi-tiNtPITONWI 
NNlPMNMI'NiOAt 
VtP«VPT*PkOAT«NlO-tl*OE 
EOPaVPT  'PA.OAT  INNI<t  >*DV 
VIP  TO  PIROT  PNtn  OP  PANIl 
EtP  10  kA«T  PNCO  OP  PANIk 
NRITI  lUPIkll  EtP>VtP«DPiNlIN 
NOITI  URPIkll  IPPlJliJaNkOiNHII 
EPT«PPT«PkOATINkIMP>t)»nP 
IP  UlTOP.OOitl  00  TO  INI 
JPNl 

00  00  J>NltNP|NOIP 

PPlJPiMPPtJI 

JP«OP*l 

00  00  JiJPiNANP 
PPlJiAl. 

JfA| 

00  too  jNHklNO.NAEf 
fPlJtlAtPUl 


JlAJtPt 

00  111  4N.)f|M0Pl 

OPIJMOi 

JElPt 

DO  too  JPNkINEOiMONEO 

PlPlJWIaPOPUl 

JOlNjPtAt 

oo  too  jijEOiMom 

ElPlJtiO. 

NkO«NlMtPTAt 

OAkk  HOONO  iTkNII  ^ 
TINPMk«TINPHk*T|N(>TIMIA 
NITUON 


PONMAT  |It|,fNlflIUt»t*tttOI 

0*0  4S 


tOhflO 

vCNNEO 

OCNNIO 

tlLMLO 

ilONEEO 

CcNNIO 

ItNETA 

tONNIO 

CeNOOO 

OOOttCC 

tOOOll 

gioDio 

II  Oil  II 
ElliONI 
IIOOII 

III  Me  01 
llOlfl 
IIOOII 
IkOOOl 
IIOlll 

loom 

loom 

loom 

llOlkl 

loom 

100  III 
loom 

IIOlll 

loom 

loom 

loom 

lOOIIO 

loom 

lOOCNI 

loom 
loom 
looiro 
loom 
Ik  01 01 
INOIII 
lOOJIO 
01)0311 

loom 

OOtlAO 

InOIIO 

lOOIOI 

000310 

OtOlOO 

OOOIIO 

lOINOO 

OOONIO 

okom 

00  ON  30 
OlOANO 
OOONIO 
OOONOO 
OOONTO 


OOONIO 

OOONOO 

OlOOlO 

ONOIIO 

OOOIIO 

IIOlll 

OIIONO 

OOllOO 

01)1100 

sooiro 

OOfllV 

OOOIIO 

looiav 

OfOOtO 

OOOIIO 

lonii 

OOOONO 

OIMIO 

IIIIAO 

OOMPO 

lOIMO 

moot 

mioe 

mrio 

OOOIIO 

oooru 

OlOTNO 

mrii 

OOITM 

ooiire 

OtOPM 

OllPfO 


J 


I ? 


u 

u 


» 

f 

M 

T 

I 


! 


3 

a 


at 

-§ 


& 

ft 


ul  g 

2!  ! 


ff  t 

Ifc  w P 


i? 

t 

t 


ss 

c 

I 


«$ 

k E 

I « 

S ' 

«• 

vt 

W M 

H 

8 


S 

B S 


-1  ^ r 

'/I  M 

teJ  K ^ 

« b »- 

>4  ui 


?»  I -j 

• ^ Id  «f 

1 i9  J r.4 

V o a »- 
O J & •«4 

O 

. J uJ  ^ O ' 
i t^  «4  I- 


• -» 2 n 

r)  u«  O U 

O ^ U i ^ ^ . 

^ Id  u;  U 

^ ^ 3 t 


>4Ulk^  *1 

lA  q U kl  3 .1)  4 It  Id  3 t 

a M ^ 5 -<  r c t-  Ub  a 1“  H 

v«e«ivs  o Mid  l»- 

MMaifftMU-n  ^Fb'yi 


i«  n:  « 

H-  ij  »»  u. 

’/  *1 

jisi^E 

M ^ J 'll 

U.  M M <« 

n H w 
U.  !.«  8*  M W9 
O iS  «r 


dlt..- 

M K M M 1 


MM  M M M M M 

X M M X M H N 

33  33333 

t«  4 J ^ 


« ^ # 4 i#  4 # 

«4  M *4  *4  «4  H H 

4 ^ ^ 4 7 

v««s|;tu^  lit 

4 # # # 4 # 7 


V > 

cc 


r » ■ » « 

hCCCt: 


!s 

«4  rt  *4  M *4  *4  M 

NS,  > N N V N 
X.  K K K K K 


«•  «l 


MiMH4$a6«:3«tbV 

is55|saf,ihss 

uC  mMCSidUtkau  out-j 


^M««*4»4«>4*4*^M«4m^*4HM  Nf4  *4 

,#«*»  9^4444444^444^^ 
I I f > I • J » • I I I I I I I J » 

»»i9xiioxinininHnm^*oin89in'*>M 

gKRRRKiitRiiSISgRiSHKIflISS; 

iSSSBS 

PwKi^^NbNNKNNNKfkKKilSNN 

«dMw4«4«4MMw4«4Mt4«4M*4M'HH«4 

aNNNsNNNN'SNVNNNsNs. 


ssiasissiassxsiisi 

PPt!£?H5sR£2tf?Sj:i!!S 
BSSSS  SSSISgBXISSeSBggSS 


i449^9t4  4 *'i4^w*4m^tf>994\ 
. ^ ^ ^ ^ ^ ^ ^ 


2« 


««  fy  H «4  ra  # 


ifl  J M I,*  k •'  MPi  f ‘ - • * 


P>  « K *4  4 *4 


S K K K Ml 


«4M^<l^^dK^«iiiyi*i«4i8n«*4MKir*4lyF^NtiiJ8tru\4u<i«|l4Mn 

asssssaaSassssassaat^^i^ssRiQiB 


iaSIliviiiiSC 

diStccf  fc  “'  “ 

N N > N Sfc  » ft  % 


41 


Table  A2.  Mwp  tor  BJRACC  Program  (CcmL  > 


-.1 


O *« 

•V  y f' 
L)  4 Ij 
U k u 
111  Ik 

k'  « 

j 

kH  vll  It 

•T  ^ 3 
*-  U 

f UJ 
Cl  Y •- 
i«  O ■« 

ty  14  V 

V a 

(/I 

> i/t^ 

•v  <y  w 

18; 

1-«  5 M 

ft  a Ki 


w i 

i Ai 

1 

1 '»•  ,.i 

u.  • s o» 

s ^ < 

M 14  • 

t U U1  / 

U k-  O 1 0 ' .1 

n Y r i-*  • 
« n (r  km  w 
(k  u.  p ri  ek  ') 
P l>  M M 

M A <y  tk  H 
a o S O'  « Ik 

.r  lb  M ,p  .1  ^4 
t-1  1/1  4 N* 

^ K a!  Cu 
.liw  u bi  m jt 

u « u > M IV  « 

jr  y u w 

^ IM  1.  «H  Lk 
D mI  UJ  ► >y  .J 
QL.  M b-  hi  W 4>  4« 
M k ^ 9t  « 4 Ui 

a ^ b*  (k  M r 

itlSS-'-'’ 

lacjkJii 

GSSSJ^'i' 


«<M««v4«4*4«4H««»4 

( < • ( I f ( < I « 
NNNNkgbJNNNN 

bifUbtblbi^btbibiM 

KSiSaKKaRag 

UOOUOOQOOO 

UOcMIUUUOUCI 

h'Kp^fktkKPk$kfkfk 

a555ia:5:;5:> 

*««4HM«4b««4e4H«4 

llSIliiSSII 

mrf/ifmwaamam 

kiK&sinls 

ttT'ji  tt'jn 


tfi  in 


W «4 

8» 


ES 


rS 

iT 

mat 


tfitfi  in  miA  tfim 


sc 


ia 

CK 


MbJMNNny  N fg  N M 


sccccs 


a 

4i 

» 

8 


KKl  K !»K  «RR!R«i^ 


If 

D O 

ou 

i4i  «4 

?- 

1/1  W1 


« 4 

. i; 

S 88 


Ittisi  !|  I 

Ul'UUi>0  uu  o 


V<  R 


y>4»  ioa>^g>Aia 

KK  »wk«V4Kk.r 

'-S^SS: 

bi  M 

'X  s 


xsxXWx  Ws  s 

mm  m 

b(  bl  ^ bi  IM  M fu  M fVi 

NVVXS*»i  XS  X 

irir  ir 

piri«nair>n  r**  m 


t 

O 

u 

14' 


R 

k 


17  w e no 
M H M M H 

cc  s sc 

in 


o o o o a q 

H M bU  H H M 

Kpeeicg 

^ VI  ^ in  M 

< I • ^ I < 


vtjii 


M in 

K e; 


4aMf4  ^ i4i0  N *^4  4«»  (4inHk4 


fj  M K »•  K v4  K 

^ *n  bj  m 


sSistiiSiKuiofiSnilHSaSnas»&Sc£lic2sKni6«lieESis&bSEE& 


41 


If  ip  for  HUtACC 


I 


•J  of  2*  ^ # 

CCSitA 

nin  nnfi 


?5?tS 

« # « •# 

»»  4 


If??? 

a**  - R 


I tsisiiiiiia 

« 

mnh 


i 


T^fale  AS.  SMopiit  Outpot  Frooi  HIRACC  Procmn  (ConL ) 


1 


\ 

t 


«0 


VO 


S 3 


M iC 
» \0 
IT 


ra 

«9 


«<a««sai<9«a9< 


— I u,  ^ . .....  . - . 


. - . ^ .)  n 

► cr  ^ 

S Ih  ^ ^ UP  « 


ti^«0^WF*U9 


4f/lUNtfS|r 

• «««««««  ««  •«•••••«•  •«••«••••• 

W M I 

r 

• -*  — - 

I I • I « I I • • I I I I I 


^ UIUJ  Ui  WIU 


I • # T I • I ♦ T ♦♦♦.♦♦♦.♦•  I 1.  » 


0\  M N i 

S wJSK  Si 

- 


JK 

« ■ 

a 


z a> 

O r* 
M ^ 

H • 

9 

s! 

S 

s 


1^ 

m 

♦ 

4 

4 

4 

4 

4 


U ca 

« ^ 

u>  Iff 

* 


i 


t 

t 


• w ^ 


K K »•  lb  K 


• #<«»«.j(y|f)^USi0K*(;ni4 


«4  N 
f4  14 

# 4 

N M 


W K 
¥r  • 


60 


AVERttE  «0TM  * .83BEB  NO.  LINES  = 12BI  T3TBL  NO.  OF  LInES=  1263 


AVRKGE  WOTH  s .03*5%  NO.  LINES  = IC*J  T3T«L  NO.  OF  LINES=  ?2«3 


I 


! I 


I 


^ j.' 

■t'j 


a ^ 

^ 4 


to 

to 


ss^ 


o 


c9o4n<Btaiait3i9«iiat«iSMi«paiVaiiB0MOriciCfr*a>w«itoo*Bj<P 


I tt  « O M 

^ 9*  9^  m ^ 


« • • • 


- - » ^ cr  T 


•JK 

U 

3 M 


m 

h 

C 


^^  F* 
(ft  itf  *0 

■ t «*  • t • 

-.  w_  - 

WfOFJWF>fOWWW#OF)>OwF>»OWWWloF»FFlW 


O ISJ 
< w> 
IJ  0) 
tf  • 


555 


Kit 

H>0 


_ . .^  (N  N N iv  M M N n FI  F 


X 

& 

s 

u 

li 

4 


5S 


1 


lUieAS.  SsoqJeOa^xitFnKn  HIBACC  nngmn  (Coot.) 


Stf^«|w90»ni«4M  ^ «0K 

• •••••»•»•••••«*•••••••  •• 

w4w4wim4 


l•«ti«Vt»«l■l«l«■I^IVC9l 


v4^Hv4«4«4|ftWv4«i4•^*4v4«4<H¥4«^^i^i^J^i^ilM^itSlM^l^J^INM^J^i^) 

iiiiir^^iiiiitfitTiTiit^i.  I iTiiii  II 

% !i;  ill  it!  iil  fe  t}!  le  II!  ^ Si  It;  iij » lil  le  % ^ lil  ^ It  iii  it*  *jj  it!  Ill » ^ 

■"■“*■■  •-  — — — • ^ ■ In  ^ lb  m •M  Ki  IC  ^ 


Wv44l^•9wi^N«4«*9N^^V«rKl^^'^0tfiONI*l0^«K 
M ( « H «4  # •«  V>  fW  K v40»  iT  IM  ^ •>  « « «!  b «0  (T  9 


5 4 S ^ ^ J J J 31 5J  5f  5, 51  ly  51 5 5. 5 3 5 ; 


kTEMCE  IIICTM  = .0385?  MC.  LIMES  = 1899  TOUL  MO.  OF  LlNES=  ‘•633 


Table  A4.  Data  Cards  for  HIRACC 


Appendix  B 

A llattng  of  « Fortran  progrant  to  plot  tho  rtauUa  from  HIBACC  In  tranimla- 
alon  U glvan  In  Appandlx  B.  Thin  program,  antltUd  TPLOT,  was  written  for  tha 
CDC  6600  and  tha  plotting  Amotlona  art  ai  Molflo  to  that  oomputar.  Thla  plot  pro- 
gram la  conalatant  with  tha  aampllng  orltarion  and  tha  tnt«rpoUllon  mathoda 
daaorlbad  In  thla  report.  Tha  four  point  Lagrango  Intarpolatlon  la  uaad  to  ganarata 
thraa  Intarmadlata  polnta  for  aaoh  output  Interval  of  HIBACC.  Tha  Hating  of  tha 
program,  tha  output,  and  tha  Input  data,  ara  glvan  In  Tabla  Bl.  Two  data  card* 
ara  raqulradi  tha  flrat  oontalna  90  eharactara  of  Idantlfloatlon,  and  tha  aacond 
ooutalna  tha  upper  nnd  lower  llmlta  (om*^)  of  tha  wava  number  range  anU  tha  dla- 
tanoa  la  Inohaa  (OX)  over  which  thla  wavenumber  range  la  to  be  plotted. 
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Tthla  Bl.  LUtlng  of  TPLOT  Program  <Coot. ) 
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